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Abstract 

The purpose of this pubhcation is to derive and discuss equations of mo- 
tion of affinely rigid (homogeneously deformable) body moving in Eu- 
clidean space of general dimension n. Our aim is to present some analyti- 
cal methods and to discuss geometric structure and invariance properties 
of the theory. Some perspectives of further developments are also dis- 
cussed. 

1 Introduction 

The theory of affinely-rigid bodies (shortly - afHne bodies) was formulated and 
treated in a systematic way by J. J. Slawianowski [ID] in terms of the 
modern differential-geometric and algebraic methods. Earlier, the concept was 
implicitly present in Eringen's theory of micromorphic continua. Those con- 
tinua, roughly speaking, consist no longer of structureless material points, but 
rather of infinitesimal affine bodies, i.e., material points with attached linear 
bases [15] [T3]. This was an obvious generalization of Cosserat's theory; the 
structure elements ("grains") not only rotate but also undergo homogeneous de- 
formations. We are aware also of the existence of some much earlier Russian 
papers (private information), however, having no explicit references we were 
unable to find them. Later on the subject developed rapidly both on the level 
of theory and applications and attracted the attention of many researchers [S] 
[7][n][^[T2][Il][5g[57][lH|[iai30j[3jy[3a[3g 

[46] [52] . What concerns theory and mathematical foundations, deep studies by 
Burov and Chevalier [S] [7] are recommended. 
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Affinely rigid body is continuous or discrete or even finite system of material 
points constrained in such a way tliat all affine relationships (but not necessar- 
ily metrical ones like distances and angles) are preserved during any admissible 
motion. Such a body performs rigid rotations, translations and homogeneous de- 
formations. It has a finite number of degrees of freedom, however, deformations 
being admitted, it is something between (metrically) rigid body (gyroscope) 
and general deformable continuum. Let us mention that some important gen- 
eral problems concerning constrained continua were studied systematically by 
Wozniak and Kleiber [H] [53] [Ml- 

This is an interesting model from the purely mathematical point of view of 
rational analytical mechanics, but there are also remarkable applications in a 
wide variety of problems. Let us quote some of them: 

• problems of macroscopic elasticity in which the length of excited waves is 
comparable with the linear size of the body 

• molecular dynamics, molecular crystals, fullerens 

• dynamics of inclusions, suspensions and bubbles 

• structured continua 

• nuclear dynamics 

• astrophysics, geophysics 

• robotics 

• modified finite elements method 

• dynamics of one-dimensional chains/lattices. 

Let us also stress some important theoretical point. Namely, in sufficiently 
small regions every smooth deformation is homogeneous; this is just the very 
idea of differential calculus. Analysis of affine bodies contributes much to the 
general theory of non-constrained continua. Any deformation gives rise to some 
field of homogeneous deformations. And fields of homogeneous deformations, if 
not necessarily holonomic, are useful when analyzing internal stresses [17j. 

Below we discuss certain dynamical models of affine bodies including prob- 
lems of the partial separability and integrability. Presented are ideas and models 
which are rather new, up to our knowledge never touched in literature. Namely, 
in all papers we know (many of them quoted in the references) it is only kine- 
matics that is based on affine geometry, but on the dynamical level affine sym- 
metry is broken and restricted to the Euclidean group of motions or to some its 
subgroups. Unlike this we discuss also models the dynamics of which is affinely- 
invariant. This is quite a new program, very interesting at least from the purely 
academic point of view. But there are also some reasons to expect that the 
suggested models are also dynamically viable and that on the fundamental level 
of physical phenomena the "large" affine symmetry of dynamical laws is more 
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justified and desirable than the restricted invariance under isometries. More- 
over, we show that instead using some exphcit potential energy expression, one 
can encode the dynamics of elastic vibrations in an appropriate from of affinely- 
invariant kinetic energy. The resulting geodetic models in a sense resemble the 
procedure of Maupertuis principle, where the orbits of motion are geodetics of 
an appropriate metric tensor build as some conformal modification of the "true" 
geometric fundamental tensor There is also some similarity to the concept 
of effective mass known from solid state physics [201 . 

In a sense of our affinely-invariant geodetic models may be interpreted as 
a discretization of the Arnold description of ideal fiuids in terms of geodetic 
Hamiltonian systems on the group of volume-preserving diffeomorphisms [5]. 
This is a very drastic discretization, reducing the continuum cardinality of de- 
grees of freedom to the finite one, namely n{n + 1), n denoting the dimension 
of the physical space (3 in realistic models) . 

It is interesting to mention about certain interesting procedures based on 
modified finite elements methods following the ideas of M. Rubin and his co- 
workers [3Sj In its most traditional version the idea of finite elements 
was to cover the bulk of deformable body with a mesh of tetrahedrons or par- 
allelepipeds (or sometimes-other simple figures) of sufficiently small size; sim- 
ilarly, the surface was replaced by a mesh of triangles or parallelepipeds (or 
other simple two-dimensional figures). Those small parts were assumed to be 
(in a good approximation) homogeneously deformable. Basing on this assump- 
tion one replaced a system of partial differential equations by a discrete system 
of difference equations. And this approximation enabled one to perform di- 
rectly the computer-aided numerical calculations. The time variable is then 
also discretized and its continuity is replaced by some mesh. But there is also 
another procedure, when instead of the spatio-temporal finite elements one uses 
also spatial ones. This hybrid method consist in an approximate representa- 
tion of continuum as a system of mutually interacting affine (homogeneously 
deformable) bodies. But the time variable is continuous and for the mentioned 
system the usual methods of analytical mechanics and qualitative theory of dy- 
namical systems are used. In various problems this hybrid discretized-analytical 
approach seems to be more appropriate. In any case the powerful methods of 
dynamical systems theory may be used. 

2 Some basic concepts. Review, notation, geom- 
etry. 

Although, as mentioned, we admit discrete or even finite systems of material 
points, it is more convenient to use the standard terms of continua. As is well- 
known, configurations of continuous media are described by diffeomorphisms $ 
of the material space N ("Lagrange variables") onto the physical space M ("Euler 
variables") [H] [TO] [TT] [H] gS] [SD] [35]. Those are assumed to be affine spaces 
of the same dimension n. Obviously, in applications n — 3, sometimes 2 or 1, 
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but for many reasons it is convenient to admit the artificially general dimension. 
What concerns fundamental concepts and notations in afhne geometry we follow 
mainly the standard treatments like [3] g], cf. also [T^ [TU] [55] [iP] [12] [15]. 
There are problems where instead of the total N we take some subset of N, 
e.g., in problems of shells, roads, membranes and strings. Sometimes N from 
the very beginning has the lower dimension than M, however, here we do not 
consider such situations; for some examples see, e.g. [32]. Motion is described 
by the time dependence of Usually (not always) M and N are endowed with 
flat metric tensors g, ij in M , N ; g e V* ^ V* , rj ^ U* (E) U* where V, U are 
respectively linear spaces of translations (free vectors) in M, N, and V*, U* are 
their duals (spaces of linear functions respectively on V ^ U). It is not necessary 
[T5] |14| . but in this paper sufficient and convenient to use rectangular 
coordinates a^, j/' respectively in N and M. They are what is often referred 
to as Lagrange and Euler coordinates. The components of 5, rj are denoted by 
ffy , VAB', obviously in rectilinear rectangular coordinates we have: 

9ij = ^ij, VAB = Sab- (1) 

The contravariant invariant metrics g^^ (1 V ®V, t]~^ € U ®U have com- 
ponents traditionally denoted by 5'^, -q"^^ (upper case indices), where 

3*V,=<5^, V^'^VCB^S^B- (2) 

Analytically $ is represented by the dependence of j/*-variables on a^- 
variables, 

y'^y\a^). (3) 

When dealing with motions, i.e., time-dependent $ (<), we introduce in ad- 
dition the time-dependence to ([3]), 

f = y\t. a''). (4) 

The placement, i.e., deformation gradient, is analytically represented by the 
matrix x with components 

^'^ ^ It. (5) 

This is a "doubled" quantity in terminology of Schouten and Veblen; for the 
fixed a G N, y {a) G Af , the tensor identification of x(a) is as follows: 

xia)eV(g)U*. (6) 

Therefore, nothing like the symmetry or antisymmetry concept does apply 
to the object %. 

Solving ([3]) or Q with respect to Lagrange variables, 

-^-a^(y*), a^ = a^(t,y'), (7) 



a 



we obtain the analytical description of the inverse mappings $ ^, ^{t) ^ of Af 
onto N. 
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Lagrange velocity field is defined on N and takes values in V, 

V{t,a) = ^it,a). (8) 
Euler velocity field is defined on M and takes values in V, 

v\t,y) = ^{t,a{t,y))=V{t,a{t,y)). (9) 

where the partial differentiation d/dt acts only on the first argument, not on 
the i- variable in a{t, y). One can easily show that: 



dy^ , ^da^(t,y) , , , ,, da^(t,y) , ^. 
v\t, y) = -g^ it, a) ^ = -XV {t, a{t, y)) (10) 



In certain formulas we need locally co-moving, i.e., material expressions of 
those (and other) quantities. They are given respectively by: 

V^(i,a) := ^{t,a)V{t,a)^x-Ht,af,V{t,a), (11) 

i^it^y) ^{t,a{t,y))v^{t,y)^x-Ht,a{t,y)f,v^{t,y). (12) 

We follow the conventions and terms used in the classical treatises by Eringen 
[14] [15], so e.g., the Green and Cauchy deformation tensors G, C are meant in 
the convention: 



Oy^ Oy-^ 

GkL = 3^JQ^Q^^9^3X'KX'L, (13) 



to avoid the crowd of characters we do not write arguments. This does not 
seem to generate confusion. Using the modern notation of differential geometry, 
namely the pull-back symbols, we write |22j : 

G = $*5, C^^-'^^'ri. (15) 

We use also the contravariant tensors G^^, G^^, where, by definition: 



Obviously 



G-i Gal = <5^l, C'^ " G,, = S). (16) 



In terms of the push-forward notation: 

G-i = $-V\ G-i = $.ry-i. (18) 
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Remark; a very important warning: There is a danger of confusion when the 
Schouten kernel-index convention is used very automatically, without sufficient 
caution. Namely, it is strictly forbidden to confuse P7|) with the metrically- 
raised contravariant versions of G, C, 

G^^ = r^i^^r^^^GAB, C'' = g'^'g^'Cai,. (19) 

One uses also the following mixed tensors: 

G^'l-.^V^'^Gal, C\ := g^'^Cai. (20) 

The tensorial "membership" of those objects (when taken at a fixed argument 
value) is as follows: 

G<eU*®U\ G-^eU®U, C(^V*®V*, C-^eV®V, (21) 

G eU ®U* L{U) , C eV ®V* L{V) . (22) 

Obviously (pij) are symmetric and positively definite. 

The pairs of tensors {G,i]), {C,g) give rise to deformation invariants. They 
may be chosen in infinity of ways, for example we can take: 

ICa M = Tr (g^) = Tr (c-^) , a = 1, . . . , n. (23) 



In the non-deformed situations we have G = j], C — g. Sometimes it is 
convenient to use the deformation measures vanishing in the non-deformed state; 
the most important of them, denoted by E, e are in ^14j referred to as Lagrange 
and Euler deformation tensors, 

E:^^iG-v), e := i (g - G) . (24) 

Obviously, their tensorial "membership" is i? G U* ® U* , e £ V* (S^ V* 
and they are symmetric. Following (|20p we can construct the mixed tensors 
EeU®U*,eeV®V*, namely 

E^B = V^'^EcB , e) = 3*'= ekj , (25) 
and the corresponding invariants 



Tr 



(^"), Trin; (26) 



obviously, those invariants are functions of (j23p . 
Spatial gradient of the Euler velocity field 
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is very important in hydrodynamics, nevertheless, it is also of some relevance for 
elasticity, and in any case it is generally well-defined in any infinite continuum. 
One can also introduce its material representation: 

Remark: as mentioned, we are using only Cartesian rectangular coordi- 
nates. Nevertheless, the above formulas remain valid in general coordinates, 
including curvilinear ones, when the usual partial derivatives are replaced by 
the covariant ones. 

The tensorial "membership" of (gT]), ^ is, 

n ev (^v* L{v) , n <=,u ®u* L{u) . (29) 

The g-skew-symmetric part of ([TT)) describes the local field of angular veloc- 
ity, 

a;^- = \ (f]^- - f}/) = \ (17^ - g,u g'' n\) , (30) 
and the (7-symmetric part refers to deformation velocity, 

d^, = \ {n^, + ^}) = \ {^', + 9,k / , (31) 

the quantity very often used in mechanics of viscous fluids. Its material repre- 
sentation, 

(32) 

i.e., strain rate, may be expressed as: 

If Gab is expressed as a function of Euler variables , then the usual time 
derivative is to be replaced by the substantial one. 

An important thing is the study of transformations and symmetries. Let 
A : M M, B : N ^ N be, respectively, diffeomorphisms of the physical 
space onto itself and of the material space onto itself. Roughly speaking, they 
describe the spatial and material transformations. In particular, spatial and 
material symmetries are described in this way. Transformations A, B act on 
configurations $ respectively through the left and right superposition: 

Ao^oB^{LAoRB){(t>)^ [Rb o La) {4>) ; (34) 

in the last two symbols La and Rb are to be meant as left and right regular 
translations by A and B. In particular, the homogeneity and isotropy (transla- 
tional and rotational invariance) of space and material are described in terms 
of such mappings. 
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3 AfRne bodies. Basic geometry, symmetries and 
canonical formalism 



Now let us go definitely to our main subject, i.e., to afEne constraints. This 
means that the set of all a priori admissible configurations $ is confined to 
AfI{N,M), the n{n+ l)-dimeusional manifold of affine isomorphisms of TV 
onto M. To be more precise. AfI{N,M) is a non-connected open subman- 
ifold of Af (TV, M) . Obviously, there are two open connected components in 
AfI(N,M). They are mutually disjoint, nevertheless they infinitesinially ap- 
proach each other along A f sing (N, M) C Af {N, M), i.e., along the closed sub- 
set of A fl {N, M) consisting of dimension-degenerating affine mappings of N 
into M, i.e., such ones which "glue" different material points. Such transforma- 
tions are non-admissible in continuum mechanics, nevertheless, at least some of 
them may be acceptable in mechanics of discrete affine bodies. Nothing bad 
happens if, e.g., all atoms of a four-atomic molecule in three-dimensional space 
happen at some moment co-planarly or collinearly placed. However below we 
shall not discuss such details, although they are interesting in themselves. If ori- 
entations of affine spaces A'', M, i.e., those of linear spaces U , V are fixed, then 
the mentioned two connected components of AfI{N, M) consist respectively 
of orientation-preserving and orientation-reversing affine isomorphisms. Only 
if we assume that N = M, then the two connected components are defined 
intrinsically, without fixing orientation. 

If we use rectilinear, in particular rectangular coordinates , respectively 
in N and M, then affine transformations, i.e., configurations of affine body, are 
analytically described in terms of linear-nonhomogeneous transformations, 

y' = x' + v'Ka^. (35) 

The quantities and ip^ k are generalized coordinates; motion is described by 
their time-dependence. 

For the sake of completeness, let us remind that <I> : A'' ^ M is, by definition, 
an affine transformations if it preserves all affine concepts, like straight-lines, 
parallelism etc. This means that there exist a linear mapping L]<^] : U ^ V , 
denoted also for obvious reasons by DL, such that 

^{a)^{h\ = ^ (36) 

for any a,b e N. Obviously, the arrow symbol denotes the vector ^ U 
originating at a S A'' and terminating at b & N. For brevity, the same symbol 
is used in (M, A'^). It is clear that for any chain of affine spaces we have 

L[$io$2] =i[$i]i[$2], (37) 

assuming of course that superpositions are well-defined. More generally, 

L[3>i o $2 o • ■ • o = L[^i\L[^2\ ■ ■ ■ L[^k] (38) 
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under the same assumption. 

We shall use the standard abbreviations GAf{M), Af{M), GL{V), L{V), 
respectively for AfI{M,M), Af{M,M), LI{V,V), L{V,V). Then, obviously, 
the operation L acts as a homomorphism of GAf{M) onto GLiV) (epimorphism 
of GAf{M) onto GL{V)). Its kernel is identical with the group of translations 
T{M) C GAf{M). Let us remind, this group is isomorphic with the linear space 

V as an additive Abelian group and consists of transformations ty : M — > M , 

V €V such that 

i„(a) = 6 if d) = v. (39) 

Obviously, T{M) is a normal divisor of GAf{M). The connected compo- 
nents of unity, i.e., subgroups of orientation-preserving isomorphisms are de- 
noted by GAf~^{M), GL^{V); their cosets orientation- inverting isomorphisms 
are denoted by GAf^{M), GL^{V) (obviously, they are not subgroups). 

Analytical representation preassumes some choice of affine reference 
frames (D; Ei, . . . , E2), (0, ei . . . e„) respectively in N and M. The points D G 
N, G M are origins and Ea G U, 6 V, A — l...,n, i = l...n are basic vectors. 
The elements of dual bases will be denoted as usual by E^ G U* , e"^ & V* . The 
parameters refer to translational motion. They are just coordinates of the 
current position $ (D) G M of a fixed material point D G N. The variables (p''K 
are generalized coordinates of relative/internal motion. When D G N is fixed 
once for all, then our configuration space AfI{N, M) may be identified with the 
Cartesian product 

Q^Mx LI{U, V) = Qtr X Q,„t (40) 

where, obviously LI{U, V) denotes the set of all linear isomorphisms of U onto 
V; it is an open submanifold in L(U,V), the linear space of linear mapping 
of U into V. Concerning the term "all" used above, again the previously men- 
tioned provisos are maintained, depending on whether the continuous or discrete 
system is meant. 

Usually D is chosen as the center of co- moving, i.e., Lagrangian (therefore 
constant) mass distribution. If /x is a positive measure on N describing this 
distribution, then O £ N is uniquely defined by 

J ol d^iiA) = 0. (41) 

Then :— $(0) G M is the current position of the center of mass in M, 
thus it satisfies: 

J o^d^^(2/) =0, (42) 

where fi^ is the 0-transport of the measure /i to M. 

Remark: It is easy to commit a dangerous mistake. Equations (HI]), (|42p 
are equivalent only for affine configurations If $ is non- affine, then in general 
the current (Euler) center of mass 0$ G M as defined by (|42|) does not coincide 
with $(i3). This is just one of definitions of affine mappings: center of mass is 
their invariant. 
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We use the (more or less) standard symbols UAf{M), SAf{M), UL{V), 
SL{V) respectively for the unimodular affine group of M, the special affine 
group of M, the unimodular group of V and the special linear group of V. Let 
us remind that (by definition) the unimodular subgroups preserve the volume, 
and their special 5'-subgroups preserve the volume and orientation. 

The fixed metric tensors rj € U* ® U* , g G V* ® V* give rise to the isom- 
etry subgroups E(N,t]) C GAf{N), 0{U,r)) C GL{U), E{M,g) C GAf(M), 
0{V,g) C GL{V). The symbol "£"' above refers to the Euclidean group, respec- 
tively in the 77- and g-sense. The corresponding matrix elements oi a — L[A], 
b = L[B], cf. (1311) satisfy then 

gij = gu a\ Oj, riAB = Vkl h^A b^B, (43) 

i.e., a e GL{V),b G GL{U) are respectively (7-orthogonal and ry-orthogonal linear 
automorphisms. As usual, the connected components of group identity, i.e., the 
subgroups of orientation-preserving mappings, are denoted by 

SE{N,r^)^E+{N,r]) , SE{M, g) = E+ {M, g), 

50({/,ry) = 0+ ((7, 77) , SO{V, g) = 0+ {V, g). 

And their cosets consisting of orientation inverting mappings are so denoted 
with the label "— " instead Let as mention, this time the and "— "- 

components are finitely-separated in the general affine and linear groups, not 
infinitesimally as previously. 

The set of affine isometrics, i.e, Euclidean mappings of (A^, r/) onto (M, g) 
will be denoted by E{N, r/; M, g). Analytically, in terms of the membership 
$ G E(N, I]; M, g) means that 

VAB = 9ij <y5 A V^B: (44) 

i.e., using the modern geometric notation, 

r] = (45) 

Obviously, becomes then 

G = 77, G = g. 

This means that is an element of 0{U,ri; V,g), i.e., it is an orthogonal 
mapping of {U, rj) onto {V, g). Just like orthogonal groups, the set 0{N, rj; M, g) 
is non-connected; it is a disjoint union of two connected manifolds. If some 
orientation standards are fixed in A^, M (in U, V), those are manifolds of 
orientation-preserving and orientation-inverting isometrics. The same concerns 
E{U, 77; V, g). The metrical rigid body configuration space is then identical with 
the manifold of orientation-preserving isometrics. 

In mechanics of affine bodies, when (f) G AfI{N,M), the transformation 
group (p4| is restricted to one with A, B being elements of GAf{M), GAf{N) 
respectively. And similarly, when some of the afore mentioned constraints are 
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imposed on then in ([M)) transformations A, B are confined to run over 
appropriate subgroups of GAf{M), GAf{N). 

We are particularly interested in geometry and dynamics of internal degrees 
of freedom, i.e., in what is going on in the manifold Qint = LI{U, V). There we 
are dealing with GL{V) and GL(U) acting as internal transformation groups. 
In analogy to ([M)) we have then the following transformations of Qmt generated 
by a e GL{V), b e GL{U): 

ip ^ aipb. (46) 

This is just exactly a consequence of ([M)) obtained by the application of the 
projection L : Af{N, M) — >■ i(/7, V) a, b being images of A, B. 

Affine constrains simplify the structure of many expressions used in mechan- 
ics of continua. The placement ([5]) becomes constant both in the physical space 
and in the body; it equals the internal part of the configuration, 

XV = V'V. (47) 

Deformation tensors are also spatially and materially constant, for example 
the expressions (ITS)) . (fT7|) for the Green and Cauchy tensors become: 

GkL = gij </5V ip^L, Cij = T]KL V^'^ ■^j, 

(48) 

The quantities fl, 57 became what in our earlier papers [18] [39] [40] [42] [43] 
was called affine velocity, respectively in spatial and co-moving representation. 

Expression ^ becomes the following affine ("linear-inhomogeneous") func- 
tion of the Euler variables y*: 

v'{t, y)-^ + ^^~' iv' - - ^(iO^ + {V' - , (49) 



il j = —j^ if j, (50) 



therefore 

compare this with (|27|) . In our earlier papers the quantity (|50p was referred to 
as affine velocity; Eringen in his papers about micromorphic continua used the 
term "gyration". Obviously, O as defined in (|50p . i.e., in coordinates- free form, 

(51, 

is an element of the tensor space V^V* ^ L{V); il^j are spatial, or laboratory- 
referred components. The co-moving (material) representation, mentioned in 
for general (not necessarily affine) bodies is given by 

n = tp-^^^ip-^nipeU(i)U* ~L{u), (52) 
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i.e., analytically, 

= <^ i = <^ i f^j '/s-'ij. (53) 

In spite we are dealing here with the special case of ([TT)) (pS)) . the quantities 
ri, are geometrically very peculiar. They are elements of the Lie algebras of 
G'L(y), GL{U). When affine constraints are replaced by stronger ones, ruled 
by some subgroups of GL{V), GL{U), then ft, il become elements of the cor- 
responding Lie subalgebras. For example, if we take Qint — 0{U,r]; V,g), or 
rather, its connected component, i.e., the configuration space of the usual rigid 
body (gyroscope) without translational degrees of freedom, then J7, f2 are re- 
spectively g-skew-symmetric and ?7-skew-symmetric angular velocities, 

= = g'" n%, h^B = -^B^ = -VBC V^"" ^""d; (54) 

they are elements of the commutator Lie algebras of SO{V,gy, SO{U,rf)' of 
SO(V,g), SO{U, Tj). Another important example is that of incompressible afHne 
body, when the internal configuration space Qint consists of volume-preserving 
(and in the continuum case orientation-preserving) linear mappings of U onto 
V . This manifold is invariant under ()46|) with a, b restricted to the subgroups 
U{V), U{U), or, in the continuum case, to SU{V)', SL{Uy. Then Vt and Vl as 
elements of Lie algebras SLiVy , SL{Uy are trace-less, 

Tr{Vt) = n\ = 0, Trin) = Vl\ = 0. (55) 

As all the mentioned groups are non-Abelian (except the planar orthogonal 
group when, n ^ 2), the objects f2, 51 are non-holonomic velocities in the sense 
of Boltzmann, i.e., they are not time derivatives of any generalized coordinates. 

In certain formulas it is convenient to use the co-moving representation of 
translational velocity, namely the quantity v{tr) G U with components 

y(tr)^ ^e^vitry ^e^^. (56) 

Canonical momenta conjugate to generalized velocities v{try and V\ = 
df'^A/dt, are denoted respectively as p{tr)i, their tensorial membership 

p{tr) € V\ P € U (g>V* ~ L{U,V), and the duality oi V G V (g) U* and 
P E U (E) V* is meant in the sense: 

{P,V) =Tr{PV) = P^,V\. (57) 

Canonical translational momenta conjugate to e J7 will be dented by 
p{tr) G U*; their components are given by 

p{tr)A^p{tr),ip'A. (58) 

The non-holonomic momenta conjugate to ft, fl will be denoted respectively 
by S G L{V) ~ ★ ~ L{V), S G L{U) ^U(g)U*. They are given by 

= P^j , t\ ^P^,^'b= ^'j f'B ■ (59) 
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Obviously, the duality between S, S and Q, fl is meant in the sense of trace: 
(S, n) = (y., = Tr (S VL) = Tr {t, ^ T.) Vl\ = S-^s VL^a- (60) 
The transformations (j46|) act on the above quantities according to the rules: 

Q.^aQ.a-'^, S^aSa^i, Q.^h-^Q.h, Y,^h-^th, 

(61) 

p{tr) p{tr) o a ^, p{tr) p{tr) o b; 

analytically the last formulas have the form 

p{try,=p{tr)j a-^\, pitrY^ ^ pMeb^'A- (62) 

Transformation rules for deformation tensor are more complicated and there 
are some subtle points about them. Namely, for the Cauchy tensor we have 

CMy C'Mfc/ a^^'^i (63) 

i.e, symbolically, 

C [aip] = if. (64) 

On the other hand, for a general b S GL{U), there is no explicit expression 
for C [tpb] as an algebraic function of C [ip]. But for isometries of {U, rj) we have 
obviously 

C [ipb] ^C[ip], beO{U,ri). (65) 

Transformation properties of the Green deformation tensor are dual to the 
above ones. So, for a general a e GL{V) there is no concise relationship between 
G \aip\ and G [ip] although for isometries of (V, g) we have obvious invariance rule 

G [aif] = G [(/?] , aeO{V, g) . (66) 

On the other hand, for a general b G GL{U), the following covariance is 
satisfied: 

G[pb]^L = G[p],,j,b^'Kb\, (67) 

i.e., symbolically, 

G [ipb] ^b*G [p] . (68) 

By their very definition, the deformation invariants (I23p are preserved only 
by isometries, 

ICp [ap^b] ^ICp[p], aeO {V, g) , beO {U, rj) , (69) 

but not under the larger subgroups of GL{V), GL{U). Concerning the last 
statement, an exception does exist, namely, the "isochoric" invariant /C„ [i^], 
which is preserved by with a G SL {V), b G SL (U), or, when admitting the 
orientation-inverting mappings, a, b being unimodular a Cz UL (V), b G UL (U). 
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The above transformation rules are fundamental for analysis of symmetries 
of equations of motion, Lagrangians, etc. 

Another very important tool is the system of basic Poisson brackets. First 
of all, let us mention that the quantities Ej, T,^b are Hamiltonian generators 
of P5|) . respectively for the a- and ^-transformations. Because of this they are 
called the canonical components of afSne spin, respectively in the spatial and co- 
moving (material) representation. Their doubled g-skew-symmetric and ry-skew- 
symmetric parts generate respectively the action of 0{V, g) and 0{U, rf) through 
. They are respectively components of the canonical spin and vorticity [T3] , 

S', - I]^-E/ = E^-g,fcS^/, (70) 
V\ = t\-tB''^t^B-yiBc^%9''^. (71) 

Remark-warning: Unlike the relationships between S and S, b are 
NOT co-moving components of spin, 

^-^^^S'j^B, (72) 

unless motion is metrically-rigid {ip ^ O {U,ri; V, g) is an isometry), and in the 
latter case V'^b are just the co-moving components of spin. 

One can also introduce the translational, i.e., orbital afhne momentum with 
respect to some fixed spatial point o G M and the total afhne momentum with 
respect to that origin. They are respectively given by 

A (o)', = x> (tr)^ , J {or, ^ A (o)^- + S^-; (73) 

when there is no danger of misunderstanding, the "label" o G Af is omitted. 

Obviously, the quantities J (o)*j are Hamiltonian generators of o-centered 
afhne mappings (ones preserving o) p4p with B — icLm analytically given by 

V = A^y'', i.e., 'x'=A',x\ '^\ = A',^k\ (74) 

obviously, it is assumed here that coordinates vanish at the fixed origin o £ M. 

Translations in M are generated by canonical translational momenta p (i?"), 
as Hamiltonian generators. Let us now quote the mentioned basic Poisson brack- 
ets 

= (5^, E V 5*^, E\ (75) 

^\,YFb\ = S^B^^D-S^D^'^B (76) 

{e^j,E^b} = (77) 

[t^B,Pc} = S^cPB (78) 
{J),Pk} = {K),pk} ^S\pj. (79) 
If F depends only on the configuration variables (a;% <pV), then 
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and if both F, G depend only on configuration, then obviously 

{F,G}=0. (81) 

Using the basic rules, the standard properties of Poisson bracket as a Lie 
bracket, and an additional important property: 

{f{K),G} = f{K){K,G}, (82) 

one can easily calculate any Poisson bracket and write down equations of motion 
as a system of equations of the form 

^ = {i^,i^}, (83) 

where H is a Hamiltonian and F runes over some maximal functionally inde- 
pendent system of phase space functions. This way of deriving equations of 
motion is much more effective than the direct use of Lagrange equations. Ob- 
viously, only non-dissipativc, variational models may be studied in terms of 
the standard Lagrange or Hamiltonian formalism. Nevertheless, once derived 
in this way, equations of motion may be easily generalized to ones admitting 
dissipation, by introducing some more or less phenomenological friction terms. 



4 AfRne dynamics. Generalized forces, balance 
laws, d'Alembert principle, afRne Euler equa- 
tions. 

Before going any further, we must introduce some additional concepts, more 
"touchable" from the practical mechanical point of view. The co-moving (La- 
grangian) mass distribution is described by the time-independent positive mea- 
sure /i on the material space N. Any configuration $ : — > M, not necessarily 
affine one, gives rise to the current (Euler) mass distribution in the physical 
space M. It is described by the positive measure /x$ on M, the ^-transport of 
/u. Obviously, it is time-dependent because $ is so. The total mass of the body 
is obviously given by 

m = iJ,{N) = J diJ,{a) = n<^{M) = J diJ,^{x), (84) 

N M 

roughly speaking, the monopolc moment of /i or Higher-order material 
(Lagrangian, co-moving) moments of inertia are given by the family of constant 
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tensors in A^, J(k) E ; analytically they are given by 

k 



J(l)^ = J a^dnia) 
J(2)^^ = fa^a^dfi{a) 



(85) 



All those tensors are synimetric. There are three iniportant special cases: 
J{o) ~ m, J(l) G U, J(2) ®U®U. Usually, and so we do in this paper, the 
origin of material (Lagrange) coordinates is chosen in the reference center of 
mass, so be definition, cf. also (PT|) 

J(l) = 0. (86) 

All the inertial tensors quoted here contain some important information: 
roughly speaking, in typical simple situations, the knowledge of the system 
of all J(fc)-s is essentially equivalent to the knowledge of fi. Nevertheless, in 
the special case of affine motion it is only J(2) £ U ®U that is relevant for 
dynamics. Then it is denoted simply as J without the label "2" and referred 
to as the "inertial tensor". It is not the same what is called "tensor of inertia" 
in mechanics of rigid bodies, however, they are essentially equivalent concepts. 
More precisely, they are linear functions of each other. I describes the inertia of 
rotational and homogeneously-deformative modes of motion. In some formulas 
we need also the covariant inverse of J denoted by J^^ eU* ®U*, 

J'''' J-^CB = 5%. (87) 

Remark; a very iniportant warning, like one concerning psp : Do not con- 
fuse J^^ with Jjj £ U* £ U* obtained form J by the 7y-lowering of indices, 

Jvab-VakVblJ^^- (88) 

In general 7^ J~^. In certain formulas one uses also the mixed tensor G 
U ®U* ~ L{U) obtained from J by the ry-lowering of the second index, 

j/b J-^^jycB- (89) 

One must be careful and avoid confusing J~^, J, Jrj, Jjj- An important 
points is that J~^, J are purely afhne concepts, whereas J^, are partially 
metrical. 
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In certain considerations one uses also the Eulerian version of J. More pre- 
cisely, there are two such versions: J[(/), o], J[(p], both being symmetric tensors, 
elements of V ®V. Obviously, o G M denotes here the fixed origin in M given 
by coordinates = . J[4>, o] denotes the Eulerian inertial tensor with respect 
to the fixed spatial origin, and J[(p] is the Eulerian inertial tensor related to the 
current position of the center of mass in the physical space, G M. Let us 
remind that in affine motion (and only then) 00 = 0(f). It is easy to show that 
the following formulas hold: 

JM'J" = ^'a^'bJ^'', (90) 
J[$,o]*^' = mx'x^ +J[<fY^. (91) 

Unlike J, which is constant and depends only on geometry of the mass 
distribution in the body, J[ip\ is configuration-dependent, therefore, also time- 
dependent. Nevertheless, J[(p] is sometimes useful as a subsystem of alternative 
generalized coordinates. 

Eulerian multipoles may be also introduced for higher values of k. However, 
for affine bodies we do not need them, and for non-afhne configurations they do 
not admit the nice orbital-internal splitting (pijl . 

In dynamics of multiparticle systems, including continua, it is also convenient 
to use multipole moments for distributions of other physical quantities like linear 
momentum and forces. When dealing with affine motion we need only monopole 
and dipole moments of those distributions in both Euler (spatial) and Lagrange 
(co-moving) representations. 

In a general non-constrained motion the l-th order multipole moment of the 
distribution of linear momentum, calculated with respect to the fixed spatial 
origin o € M is analytically given by 

JCoiiy'-"' = j y'^ ...y'^vHy)dii^{y)^ (92) 
= Jy(ay'...y{ay'V^a)dfi(a). 

As usual, denotes the 0-transport of the measure /i from N to M, and v, V 
are, respectively, the Euler and Lagrange velocity fields ^ ([5]); the coordinates 
are assumed to vanish at the origin o G M. The tensor JCo{l) is symmetric 
in the first /-tuple of indices. It depends explicitly on the choice of o G M. One 
uses also quantities JCinti}) for which the fixed reference point o is replaced by 
the current position of the center of mass o^,, 

/c™t(0'"-"' ' = J {y'' ~ ■ • • iv" - ^'iy)df^4y)- (93) 

Being tensor in V, JCo{l), /Cmt(0 Euler-like quantities. In certain formu- 
las it is more convenient to use the multipole moments with respect to Lagrange 
variables a^; we denote them by IC{1). 

jC(l)Ai-Au ^ f a^' ...a^^V'{a)dfi{a). (94) 
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This is a mixed tensorial quantity, namely, the l-th order symmetric con- 
travariant tensor in U and the usual vector in V. In spite of its being injected 
in the "abstract" material space, in many formulas IC{1) is more convenient and 
effective than the purely spatial IC{1). 

Obviously, the monopole moments are identical with the total linear mo- 
mentum in its kinematical version "inertia x velocity; we shall use the symbol: 

F := ICoiOy = ICiOy = m— = mv{try. (95) 

As is well-known, in analytical mechanics there is a subtle distinction be- 
tween canonical linear momentum p{tr)i and kinematical linear momentum fc*. 
The relationships between them is a dynamical concept, not kinematical one. 
Below we shall return to this problem. 

The dipole moments /Co(l)- l^int{i), ^(1) have to do with velocities of 
rotational and homogeneously deformative motion, combined multiplicatively 
with appropriate internal objects. The doubled skew-symmetric parts of /Co(l), 
/Ci„((l) represent respectively the total angular momentum with respect to the 
fixed origin o € M and the spin angular momentum (one with respect to the cen- 
ter of mass) , both in kinematical versions (vector product of the radius vectors 
and kinematical linear momenta). 

In affine motion we obtain for ^ = 1: 

ICJ' = m x'v{try + i^\'^J^B^ (gg) 

This is just the additive splitting into translational ("orbital") and internal 
("affine spin") parts, 

/C„t/^" = mx'^;(^r)^ /C,„/'^ = ^\^J^^. (97) 

The internal part does not depend on the choice of o. To avoid the crowd of 
symbols we denote it simply by /C*^ without the label "int". The corresponding 
translational and internal angular momenta, both meant in the kinetic sense, 
are denoted respectively by: 

^ IC.tr'' - ICotr , S'^ = /C - /C J\ (98) 

The total angular momentum with respect to o € M is denoted by 

Zo'' ^ Co'' + S'' . (99) 

Obviously, in the physical dimension n = 3, the skew-symmetric tensors (j98p 
l|99p are identified with the corresponding axial vectors with components Cg, 

The nice and intuitive translational-internal splitting like holds also for 
non-afhne bodies, but only for the dipole moments, i.e., for /Cc,(l), no longer for 
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ICo{£) with £ > 1. This fact has some deep geometric reasons. Unhke this, ah 
Lagrange moments IC{£) spht in this way. 

It is convenient to express also the distribution of forces in terms of its muhi- 
pole moments. These moments play the role of generalized forces responsible for 
the dynamics of some collective models of motion. As expected, the monopole 
and dipole moments of forces are sufficient for describing affine motion. 

Multipole moments No{l), Nint{i), of the distribution of forces are 

given by formulas obtained form ([M)) ([M|) by substituting instead veloc- 
ity field the density of forces per unit mass. This density at a given particle 

d<s> - 

dt 



a G iV and the time instant t G R will be denoted by a;$, analyt- 



ically T'^[t,a;y^,^]. The dependence on configuration $ (analytically as 

functions of a*^) and generalized velocity (^^ as functions of a'^) may be func- 
tional, in particular non-local in space and time (memory); this is the reason 
of using the brackets symbols []. However, in applications we have in mind, 
mainly simple bodies (usually elastic or viscoelastic), J^* will be of the form 

J-^ (^t, a]y^t, a),x-'Kit,a), ^^{t,a)j , i.e., local and memory-free. Quite inde- 
pendently on any particular model, the meaning of J^* is that the quantity 



dn{a) (100) 



is the force acting on the material mass element dfi{a) at the time instant t, when 
configuration is given by $. Substituting J'^ instead in formulas (|92|) (|93p 
((94|) . we obtain, as mentioned, generalized forces No{l), Nint{l)i -^(0 controlling 
the dynamics of some collective modes. And again for / = 1 (and only then) we 
have the orbital-internal splitting, 

NJ' = No tr' + N,^' = x'F^ + Nni' , (101) 

where the monopole 

^ f Pdfi{a) (102) 



is the total force acting on the system (in a sense on its center of mass). 

The quantity N^'' will be also referred to as an affine moment of forces or 
"affine torque" with respect to the origin o S M, and A^*^^ is an "internal affine 
torque"; it is related to the center of mass as the "origin of lever". Their doubled 
skew-symmetric parts are usual torques (moment of forces), 

- No'' , N,J' - N^J' ■ (103) 

In affine motion, when configurations are given by p5|) and placements by 

all these generalized forces are functions of a;*, ^; (p'^x, ^^^'^dt^ ^^'^ possibly 
of time t explicitly. 
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The power V of distribution of forces is given by 



V 



dt 



(i, a)d^la). 



(105) 



It is easy to show that in affine motion this becomes 

V = + g„,, iV™ n\ = Gab + Gcb N^"" ^""a 



where for simpUcity the arguments of F, N are omitted (cf (|104p '). 
abbreviation for v{tr)\ and denotes "gyration", i.e., affine velocity, cf. 
(|50)) . ((5T|) . For simpHcity the label "ir" at iV*-'' is omitted. 
Non-constrained equations of motion have the form: 



(106) 
is an 

(BSl), 



and the non-constrained kinetic energy is given by 

1 r dy\t,a)dv=[t,a) 



T : 



9i3 



dt 



dt 



d^ia). 



(107) 



(108) 



When affine constraints are imposed, then ?/* must be expressed in terms of 
generalized coordinates {x^.lp'^a) like in ([35|) . and, according to the d'Alembert 
principle, the right-hand side of (|107p must be modified by introducing the 
distribution J-p; of reaction forces maintaining the constraints. 



{t, a) = r 



Fr 



dy-' 

t,a; y^{-,-), -gfi-r) 



(109) 



Substituting to p08l) the parametric description of constraints, 
obtains: 



T = Ttr 



_ m dx^dx^ I dip\ dip^s jAb 

Unt- 2 9^3 + 2 dt dt ' 



(110) 



where, obviously, Ttr,Ti„t denote respectively the kinetic energies of transla- 
tional and internal motion. This splitting, i.e., the absence of translational- 
internal interference terms is based on the assumption (|86p which tells us that 
the origin of Lagrange coordinates is chosen at the reference center of mass. Ex- 
pressions (j35p . as any constraints, may be, or rather should, be automatically 
substituted to the kinetic energy (|108p . however, it would be absolutely wrong 
to substitute them to reactions-free equation of the unconstrained motion (|107p . 
They are to be substituted to ()109p together with constitutive laws for reactions. 

Constraints are assumed to be ideal (passive), i.e., the power of reactions 
Pb( (|105P with J'^ji substituted instead J^*) does vanish on any virtual velocity 
compatible with (pSj) . 



^ ^ dt ^ dt ' 



(111) 
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thus 

Vr - .9,, Fr'v' + g^j 7V«™ n\ = Gab Fr^ v"" + Gcb Nr^"" ^""a, (112) 

for any {v^,fl^j). But this means that the monopole and dipole moments of 
P09p are free of reactions; the total force of reactions and their afRne torque do 
vanish, 

Fd=0, Nd^^O. (113) 

And finaUy we obtain the effective system of ordinary second-order differ- 
ential equations imposed on the time dependence of generalized coordinates 



or, in the balance form: 

^ = ]V« + J-"" = !i',„ Jl^l'- + N". (117) 

As previously, written without any label is an abbreviation for Kinti^Y'' , 
the internal afhne momentum ("hyperspin"); similarly, iV*-' denotes the afhne 
torque ("hyperforce", "dipole of forces") Ninti^Y'' the "lever arm" of which orig- 
inates at the center of mass. 

There is an important difference between (|116p and (|117p . Namely, if the 
total force vanishes, — 0, then (|116p becomes the conservation law of the 
linear momentum (in kinematical form). Unlike this, (|117p does not become the 
conservation law for afhne spin even if the affine torque vanishes. An obstacle 
is given by the second term on the right-hand side of (|117p . If one takes into 
account (|116p . then for any choice of the spatial origin o G M, pizp may be 
written as: 

— — ^ Na + m v{tr) v(tr)' ^ —J . (118) 

at at at 

Equations pi7p (jllSp may be concisely written in the following suggestive form, 
a bit symbolic one: 

^^N'J +2^^, (119) 
dt dgij 

^^n;j+2^. (120) 
dt agij 

The parametric dependence of Tint, T on the metric tensor g V* ^ V* 
and the resulting terms on the right hand sides, just visualize the fact that the 
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non-conservation of K^^ , Kg-' , even in the interaction- free case, follows from 
the metrical breaking of afHne symmetry of degrees of freedom. And it is very 
intuitive that the internal and total angular momenta 5*-', Zo''\ i-e., the doubled 
skew-symmetric parts of K^^ , Kg-' , are conserved quantities when the corre- 
sponding affine torques iV*-' , N^^ do vanish; moreover, they are conserved when 
affine torques are symmetric tensors, i.e., when the usual torques (jl03p do van- 
ish. This is just the conservation of kinematical angular momentum, or its 
balance, when A/"'-^, AfJ"^ are non- vanishing. There is an obvious similarity to 
the relationship between conservation of angular momentum and symmetry of 
the Cauchy stress tensor in continuum mechanics. This is quite natural because 
the volume density of contact forces in continuum is given by the divergence of 
the Cauchy stress field cr*-' : 

r-crHj. (121) 
The corresponding expression in terms of Lagrange variables has the form: 

p ^T''\k^^T^\ (122) 

where T denotes the first Piola-Kirchhoff stress tensor. Substituting this to the 
definitions of N^^ , iV^' performing integration by parts and taking into account 
the conditions at infinity, one obtains: 

N^' = - J N^' = - J (123) 

f^AB _ rpAB 



where, is the second Piola-KirchhofF stress tensor. Obviously, integration 
is performed respectively with respect to the spatial and material rectangular 
coordinates. Expression in terms of general coordinates is obvious, for example 
in (|122p (I123P one must then replace the usual partial differentiation by the 
covariant one. 

It is interesting to rewrite the balance equations (|116|) (|117p in terms of the 
co-moving and mixed representations. Using the co-moving components: 



we can rewrite (|116p (|117p as follows: 

iB 7-1 T?CA I 9A T>AB T-1 T>DB , ^TjAB ^ir,r\ 

^ ~k J BcK- +i' , — — = J cdK +N . (125) 

Let us also mention the mixed, spatial-material form of the internal balance, 
i.e., second subsystems of (|117p (|125p . namely: 

dll^d^jBA^ ^M^ (^26) 

dt dt^ ^ ' 



(124) 
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In a sense, (|125l) are "affine Euler equations". As said, "in a sense" only, 
because the non-dynamical (forces-independent) terms on the right-hand side 
do not vanish even in the case of highest inertial symmetry, when 

jAB^j^AB^ (127) 

Using the co-moving components of velocities, we can write (|125p as follows: 

™ '^^^ ™ 0-4 -B I 9a dfl^c jCA oS O^' tCA I CtAB /tiq\ 

m—T— = —m\l bv + r , — - — J — —il d il c J +^ ■ (128) 

This follows from the obvious fact that the second of equations ([77)) may be 
alternatively written as: 

k'^B ^ ^B^ jCA^ J^^J ^ J[(p]™, (129) 

J[ip] given by ^U\\ . Obviously, the first formula in p29p is more convenient, 
because the relationship between and K is based on constant coefficients. 

The balance formulation of equations of motion, (|116p (|117p pisp pi9p 
(|125p is very convenient when discussing some additional constraints imposed 
on affine motion. The point is that often one means constraints expressed 
mathematically in some geometric terms. The most natural of them are: 

1. Metrically-rigid body, i.e., gyroscopic motion. Then the affine velocity fl 
is g-skew-isometric, and, equivalently, Q is ry-skew-isometric: 

^Tj+n/ = n'j+gjk^'^mg"'' = 0, (130) 

n^B + ^B-^ = h%+VBCV^''^^D=0, (131) 

i.e., O, are elements of the Lie algebras SO{V,gy, SO{U,riy of the 
corresponding orthogonal groups. And, obviously, if permanently remains 
within the manifold of linear isometrics 0{U, rj] V, g); more precisely within 
one of its connected components. There are ^n{n — 1) degrees of freedom 
of internal (relative) motion . 

2. Shape-preserving motion, i.e., superposition of gyroscopic and dilatational 
motion. This means that ip permanently remains within the manifold 
M.'^0{U,ri;V, g), of linear-conformal mappings. Therefore, il, are re- 
spectively elements of the Lie algebras R(SSO{V,gy, R®SO{U, r/Y of the 
linear-conformal groups ]R+S'0(V, g), R+S'0(J7, 77). This means that they 
split uniquely into sums of skew-symmetric and identity transformations, 

Q = uj + aldv , n = Lj + aldu. (132) 

Here w, u) are respectively g-skew-symmetric and 77-skew-symmetric, cf. 
(|130p . p3ip . and Idy, Idu are identity mappings in V, U. There are 
in(n — 1) -|- 1 internal degrees of freedom. 



23 



3. Incompressible afRne body. Then tp preserves all volume standards, and 
affine velocities SI, are traceless mappings, i.e., satisfy two equivalent 
conditions: 

n\=Q , n\ = 0. (133) 

This means that they are respectively elements of the Lie algebras SL{Vy , 
SLiJjy of special linear groups SL{V)' , SL{U). There are {it? — l) inter- 
nal degrees of freedom. 

4. Purely dilatational body. There is only one internal degree of freedom. 
The configuration space may be represented as a manifold of linear iso- 
morphisms = Xipo , where A runs over K"*" , and Lpo is some fixed isometry. 
The particular choice of ipo is non-essential. Obviously, affine velocities 
are then one-dimensional objects proportional to the identity mappings: 

n^alAv ^'^Idv , n = aldu ^'^Idu- (134) 

5. Rotation-less motion, i.e., purely deformative motion. This must be some- 
thing completely opposite, complementary, to the gyroscopic motion, i.e., 
to the item (1). The simplest possibility is to replace ()130p by the condi- 
tion that n is (^-symmetric, 

n'j ~ = n'j - g^k g"" = o. (135) 

In any case this is the most natural and geometrically unique possibility 
of defining the rotation-free motion in V. But some important novelty 
appears now. Namely, constraints (|135|) are non-holonomic, unlike the all 
formerly quoted. The geometric reason for that is such that ij-symmetric 
operators do not form a Lie algebra. Commutators of such operators are 
just g-skew-symmetric. There are no restrictions on configurations, but 
only ones imposed on virtual velocities. The arena of motion is given 
by some submanifold in the 2n^-dimensional space of Newtonian states 
which are parametrized by {ip'^A,^\) or equivalently by . The 

mentioned submanifold is parametrized by ip'^A and, let us say, independent 
components of (Sl j + f^/), i < j. This gives together + {n + 1) = 
in (3n -|- 1) independent state variables. Such a rotation-free motion may 
occur, e.g., when one deals with suspensions in a very viscous fluids. In 
analogy to the item (1) one can also try to define rotation- free motion in 
co-moving terms, i.e., as the 77-symmetry of Q, 

h\ - f^B^ = h% - l^BC ^^D v""^ - 0. (136) 

But now a new surprise appears, namely equations (I135P (I136P are non- 
equivalent, unlike (|130p p3ip . Indeed, (|135p implies the G-symmetry of 

h, 

n^B - Gbc ^^d ""-^ = (137) 
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rather then its 77-syninietry (|136p . And conversely, (|136p impHes the C- 
symmetry of O, not the g-symmetry given by (jl35p . 

r!'j-C,fer!*^™C-i™ = 0. (138) 

And again, (|138p is something else then (|130p . Both models are interesting 
from the formal point of view of analytical mechanics. But it seems that 
it is rather p35p that is physically more applicable. 

Equations of affine motion with the above all constraints 1-5 may be de- 
rived from the d'Alembert principle, which tells us that for ideal constraints 
the moments of reactions do not do any work on virtual velocities, i.e., on 
VL-s satisfying (|130p - (|138p . For gyroscopic constraints the virtual affine veloci- 
ties f2 are gi-skew-symmetric, thus the reaction moments Nr are symmetric and 
the effective, reactions-free equations of internal motion are given by the skew- 
symmetric part of the balance laws (jll7p . i.e., by the skew-symmetric part of 
(|115p . therefore, 

^-^ = iV^^-iV^^ i.e., ^=AA^^. (139) 
dt dt dt 

This is the balance law for spin. Explicitly we have 

d'^ip^B jAB i d'^f'B jAB Kfi] /i.n^ 

^■^^^ -^""^Z = ' ^^^^^ 

gtjifi^AV^B = VAB, (141) 

where the last formula is the final, integrated form of ()130p . Equations of 
gyroscopic motion consist of the joint system (|140p fc (|14ip . Let us stress, it 
would be wrong to substitute constraint equations ()14ip simply to ()115p or 

(dUD. 

The implicit equations (|141l) may be transformed to the parametric form 
where tp^A are expressed as functions of some generalized coordinates g" on 
O {u, 77; V, g), e.g., Euler angles, canonical coordinates of the first kind (rotation 
pseudo- vector if n = 3, or rotation bivector for the general n), etc. The functions 
V^A (9") niay be simply substituted to (|140p . resulting in a system of (n — 1) 
equations for in(n— 1) coordinates q" {t). Let us stress again that such a 
substitution to (|115l) (I117P would be wrong. This pattern is to be followed in 
discussion of all other mentioned holonomic constraints. 

Equivalently, the gyroscopic balance equations (|139p may be written in co- 
moving terms, as a skew-symmetric part of p25p . p28p . e.g., 

C jCA _ ^^jCB ^ _^B^ ^D^ jCA ^ ^A^ ^D^ jCB _^ (^43) 



dt dt 



f^AB _ fjBA_ 



Equations of motion of the body with frozen shape consist of (|139p / p40p 
and the g-trace of ([TT7|) / pT5)l . 

9^^^-9.N^+g,,^^J , (143) 
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I.e., 

5.,^V^ J^^=g.,iV*^". (144) 

The system (fT39)) ([T43)) . i.e., plO)) (fT44| contains in (n - 1) + 1 independent 
equations and it must be completed, e.g., by the parametric representations of 
satisfying: 

^\^\^\ , AeK+ , *eO(f/,?7;y,g), (145) 

therefore. 

Gab = 9ij 'p'a^b = VAB- (146) 
Here A is dilatational generahzed coordinate, the remaining ones may be chosen, 
e.g., as Euler angels, rotation bivector (axial vector when n = 1), etc, parame- 
terizing 4' G O (L/, ?7; V, g). This parametrization, Lp = (?"), a = (n — 1) 
may be directly substituted to the system (|139p / (|143p (|140p / (I144p quite auto- 
matically, the reaction forces do not occur there. 

Let us observe some important novelty, in co-moving representation the con- 
traction in the trace of p25p in principle is not performed with the help of the 
metric rjAB but instead, with the help of Green deformation tensor G. However, 
this does not influence anything, because in view of (|146|) one con simply divide 
both sides by and write: 

VAB ^ = -if CD i?^^ VAB + VAB iV^^- (147) 

at 

The effective, i.e., reactions-free equations of isochoric (incompressible) mo- 
tion ([3]) are given by the g-trace-less part of 11171 [TTSl 

^ (^if'-' -igafcif^V^) = (a^'^ - ^ .9at iV^V^J (148) 



B „ij \ jAB 



\ dt dt n ^""^ dt dt ^ ) 

There are (^-n? — l) independent equations in this system, and it must be com- 
pleted by substituting constraints equation: 



And again, when using the co-moving representation (I125P we must take the 
traceless part in the sense of the Green deformation tensor G, not in the sense 
of ?7, i.e., 

= (if CD if - ^ GklK^^ CD K^" ^""^ + 

+ N^B _Iq^^^CD ^-lAB 
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Just this form is compatible with the d'Alembert principle, and now the use 
of G instead rj is essential, not cosmetic. 

For the purely dilatational motion Q d'Alembert principle implies elimina- 
tion of reactions by taking as equation of motion the g-trace of (I117P / (|115p , 



I.e., 



Substituting: 



ip\ = A^-'a, Gab = P VAB = >^ Vab, (153) 
where \l/ is a fixed reference isometry, we obtain: 

?7i^L J^^A ^ = 5., iV*^ (154) 

where, obviously, after the substitution of constraints gij N^^ becomes the func- 
tion of ^p, equivalently, of (A, As usual, the explicit dependence on 
time t is also admitted. 

As usual, when we use the co- moving representation (|125p . the trace scalar 
is in principle meant in the sense of G, 

Gab-^ = -GabK^"" CD K^"" + GabN^"", (155) 
at 

but (|153p enables one to contract (I125P simply with the help of i], 

VAB^^ = -VabK^"" J-^cdK^'' + i1abN^''. (156) 
at 

Finally, for the rotation-free motion, i.e., for non-holonomic constraints (jl35p . 
d'Alembert principle implies the symmetric part of pi7p / (|115p : 

dt dt dt dt ' ^ ' 

^a^J^^ + ^^a'^ = N^^+N^\ (158) 

This is to be completed by equations of non-holonomic constraints (|135p ; reac- 
tions are automatically eliminated by the symmetrization in (|158p / (|157p . 

The purely dynamical term on the right hand sides of (|158p (|157p , N^^ + N^^ 
is symmetric. The symmetric part of the affine moment of reactions vanishes, 
therefore, reaction forces do not enter those equations, A'^ '"' + Nfj -'^ — 0. The 
skew-symmetry of Ar is due to the fact that it must be g-dual to all virtual 
velocities compatible with constraints (|135p therefore, 

iVfl'^w,, = 0, (159) 



27 



for any twice covariant symmetric tensor — ujji . 

There is however some dehcate point when transforming (|157|) (|158p to the 
co-moving representation. Namely, the twice contravariant symmetric tensor 
7\ry' -|- TVJ' is then transformed into symmetric tensor N^^ -\~ N^^, and the 
corresponding co-moving reaction moment Nn is skewsymmetric just hke the 
spatial one, Nfj^^ + NfP^ = 0. And one might suspect some mistake or misun- 
derstanding because one cannot prove that is dual to 77-symmetric co-moving 
affine velocities. But everything is correct, one cannot prove, because it is not 
true. Namely, if deformation occurs, then equation (|159l) . more precisely the 
system: 

Nd"" {gkj m) = 0, gkj n\ ~ 0.^ = (160) 
does not imply that 

Nr^^ Ucd O^b) = 0, VCD ?l'^B - VBD n^'c - 0. (161) 



The point is that the ^-symmetry of f2 does not imply the 77-symmetry of il; 
instead, 51 is G-symmetric and in virtue of (jl60p we have 

AC fn^^ f>D\ _ n n r,D^ „^^r,D 



Nr""^ [Gcd b j = 0, GcD B - VBD ^c = 0. (162) 

instead the wrong formula (|16ip . 

A similar problem appears when materially-non-rotational constraints (|136p 
are discussed. And something rather strange is obtained from the d'Alembert 
principle. Namely, the afhne moment of reactions, Nr/' must be such that 

Nd"" Qkj m = 0, (163) 

for any affine velocity J7 such that ft is ?7-symmetric, i.e., satisfies (|137p . But 
this means that ft is C-symmetric, i.e., satisfies (|138p . Let us introduce the 
mixed tensor T) & V* ® V given analytically by 

T)^ :=gk,nC-'"'^. (164) 

The duality between reaction moments Nr^ and afhne velocities fl satisfying 
(|138p implies that the tensor N defined by 

Nd^ :=Nii''Dk', (165) 

must be skew-symmetric: 

Nd^ + Nd' = 0. (166) 

Therefore, the procedure for obtaining reactions-free equations of motion is 
as follows: Take pi7p / (|115p and transform it by multiplying on the right by J): 

at at at 

J^"" ^rn' - N^'. (168) 
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Finally we take the symmetric parts of these tensor equations. In this way one 
obtains reactions- free equations of motion with non-holonomic constraints (jl36[) ; 
everything based on the assumption of validity of the d'Alembert principle. 
Physically rather strange and mathematically rather obscure model. But it 
should be investigated if the analysis is to be complete. Similarly, when using 
the co-moving representation (I125p . one should multiply it on the right by the 
CO- moving representation of (|164l) , 

:= ""j ^k' = Gac (169) 

Then one obtains: 

^ Sc^ = -X^*^ J- Vjvi^'^'^ ^c^ + iV^^ Sc''. (170) 
dt 

Taking the symmetric part of (|170p we obtain the effective reaction-free 
equations of materially rotationless motion, because the d'Alembert principle 
implies that 

iV/'^ Sc^ Nj^'^ Sc^ = 0. (171) 

5 The link between kinematical and canonical 
concepts. Hamiltonian and dissipative models. 

Unfortunately, kinematical and canonical, i.e., Hamiltonian, concepts are often 
confused in mechanics and no sufficient attention is paid both to distinctions 
and links between them. The finite-dimensional model we are dealing with here 
enables one to understand those details correctly and sheds some light on the 
more difficult problems one is faced with in non-constrained continuum mechan- 
ics. In section 3 after the basic discussion of kinematical quantities like various 
measures of deformation and velocities we introduce phase-space concepts like 
canonical translational momentum, affine momentum, affine spin, angular mo- 
mentum and metrical (usual) spin. The geometric meaning of canonical linear 
momenta, canonical affine momenta, etc. was that of Hamiltonian generators 
of natural groups of affine transformations acting in the physical space and in 
the material space (in the body itself). The basic Poisson brackets were given 
and it was seen they were expressed by the known structure constants of affine 
and linear groups. And then equations of motion of conservative systems were 
given in the canonical form (|83|) based on Poisson brackets, does not matter 
where the Hamiltonian iJ was taken from. And, as a rule, equations of internal 
motion are given by the balance law of S or E, 

= |E-e = {£Vi^}. (172) 

This is the way of thinking of theoretical physicists working in fundamental 
problems. But at the same time, it is clear that usually in mechanics of continua 
one proceeds in another way, using kinematical quantities like linear momentum 
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fc' or internal angular momentum 5*, where, according to the school wisdom, for 
the material point k'' = mv'', denoting the translational velocity, and angular 
momentum is the vector product of the radius vector and linear momentum. 
And one begins from the system of Newton equations, and constraints are taken 
into account on the basis of d'Alembert principle. The relationship between 

kinematical quantities A;'^, X*-' , X'^-^^ and canonical ones {pi,pA,'^^j,'^'^B 
is based on Legendre transformation. Denoting the background Lagrangian of 
the non-dissipative mechanics by L, we have 

Qyt QyA Q^]. dn^A A 

depending on if Lagrangian is expressed respectively on velocity arguments in 
the version w', v"^, Q^a', V^a = '■P^A- If Lagrangian has the potential form: 

L^T-V ^Ttr + T,,,t-V{x\^\), (174) 

T, Tint: Ttr givcu by the usual formula (jllOp . then Legendre transformation 
expresses canonical momenta pi, P\ as the following functions of generalized 
velocities: 

Pz^m g,j — ^ g^jk\P ^ = g.^ J . (175) 

Kinetic energy (jllOp and Legendre transformation (|170l) (|175p may be writ- 
ten in some alternative forms, very useful in theoretical analysis, e.g., 

Ttr = ""-^Gabv^v'', (176) 

T^nt = \gab?i\^''lJ'''' = \gz,n\mj[^f\ (177) 

PA - mGABv'' = GABk'', (178) 
S-^B = Gbi3 17^cJ^^, = ffjfc r!"™ JM"^ (179) 

Let us notice that the second equation of ([57)) may be alternatively written 
as follows: 

K'^ ^Vl\J[ipf\ K^^^n^cJ^'^. (180) 

therefore, 

= K^'ghj, t^B = K^'^GcB. (181) 

Now we have the complete "dictionary" between two versions of concepts 
like linear momentum, afhne momentum, angular momentum, afHne spin and 
just spin. Those are kinematical and canonical (Hamiltonian) versions. The 
relationship between them depends on the particular choice of dynamical vari- 
ational model, i.e., on the choice of Lagrangian. The above formulas are valid 
for the classical potential class of models. If we admitted in L some velocity- 
dependent terms in addition to the kinetic energy, the relationship would be 
different. For example, if magnetic forces are present, Lagrangian differs from 
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I|17ip by terms linear in generalized velocities. Then Legendre transformation 
expresses kinematical quantities as affine, i.e., linear- nonhomogeneous, functions 
of kinematical ones. Inverting Legendre transformation, i.e., expressing kine- 
matical quantities as functions of canonical ones, and substituting them to the 
energy function, 

E = T + V ^Ttr + T,^t + V{x,^), (182) 
one obtains Hamiltonian 

H = %+V ^%r+%nt+V{x,^). (183) 

Obviously, the kinetic term of the Hamiltonian is given by: 
In analogy to p76p . (|177p we can also use the following suggestive expressions: 

%r = =^G-^^''pAPB, (185) 

Am 

%nt = ij^^S^^S^iG-i^^ = (186) 

Let us observe that the quadratic forms piOp (|184p have constant coeffi- 
cients. Unlike this, ()176p (|177p (|185p (|186p have configuration-dependent coeffi- 
cients, however, the geometric objects $7, Vt, E, S are more suggestive than ^^^j^, 
P^i. The reason is that they are Lie-algebraic objects relevant for the structure 
of our configuration space; for example E j, T,^b are Hamiltonian generators of 
GL{V)^ GLiU) acting on the internal configuration space. 

Roughly speaking, equation (|18ip tells us that the kinematical and canonical 
affine spin, K^V®V,T,&V®V*-^ L{V) are related to each other by 
the gr-shift of the second index. The metric tensor g is fixed, constant and 
in appropriate coordinates its components are given by the Kronecker symbol, 
therefore, analytically this is a rather cosmetic difference. The co-moving objects 
S S U ® U* ~ L{U), K ^ U ® U are also related to each other by some 
shifting of the second index, however the shifting tensor depends on generalized 
coordinates, it is simply the Green deformation tensor. 

The most important geometric and physical content of equations of mo- 
tion is summarized in balance equations for Hamiltonian generators (p^, E j), or 

dpA 



dt 



^{PA.H} , ^ = {E-4^,i/}. (188) 



The right-hand sides may be calculated with the help of the basic Poisson 
brackets (|75 |) -([8T |l . The fundamental properties of this operation must be used. 
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e.g., Lie-algebraic rules, and first of all, the separation rule Substituting 
to (|187[) (|188[) the Legendre transformation, e.g, (|175p for the potential systems, 
or better its equivalent forms (|180p psip . one obtains second-order differential 
equations for the time-dependence of generalized coordinates (a;*, 1^9*^1) . The 
balance laws (|187p for the potential models (I17ip have the form: 



where 



NW] = {^'j,V] = -^A-^, = -9''^. (190) 



The differential operator which acts on V on the right-hand side of (|190p 
equals the minus generator of left G'L(T^)-mappings acting on the argument of 
V. Calculating the Poisson brackets |Sj-,T}, g-raising the covariant indices, 
and substituting ([T75t . (fTM|) to HHH), one obtains just (ITTB)) (fTT7|) or 

equivalently (|125p with 

AT^J' = ivyy] g'^J, TV-^^ = v'-i^i (^-i^jiV^J. (191) 

The distinction between tensors iVj [y], and N^^ is, obviously, just as one 
between E*,- and K"^^ , of a rather "cosmetic" character. Indeed, the shift of indices 
is performed with the use of a fixed tensor g, the matrix of which in appropriate 
coordinates coincides with the Kronecker symbol. (According to the commonly 
accepted convention we could use in principle the same kernel symbol and write 
simply S*^ instead K^^; we do not do it, because certain misunderstandings 
would be possible). And again some delicate problems appear when we compare 
the spatial and material descriptions. Namely, in the co-moving representation, 
(fTSg]) becomes 



dpA dV j 



dt dx^ dt 

where 



{^''B,r] + NMV] (192) 



N^^[V]^[t^s,v}^~^^B-^. (193) 

Substituting here the Legendre transformation and (|180p (|18ip . we obtain 
just the representation p25p . It must be stressed however that N^^ in ()125l) 
is not the 77-raised version of A^'^_B[y] in ()193p . but instead it is its G-raised 
version: 

fjAB ^ nA^[v] G-^^^ ^ N^c[V] if^. (194) 

And here really the use of two different kernel symbols would be more ade- 
quate, however, we are afraid of the crowd of symbols and of changing them in 
the course of writing. To obtain (jl25p one must perform some calculations, for 
example show that 

^ = J-\dK^^Gcb + J-'bdK^''Gca= (195) 
dt 

= GbC ^'"'a + GaC ^'"b ~ i^ij — ^ji) 'P^A '^"'b- 
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where, obviously, 



In analogy to this we have 



^^-(nAB + nBA)^-'^^-'''j, hAB:=VAC^%- (196) 

The proof is easy and we do not quote it here. The formulas (|195p (|196p are 
very suggestive and reveal some more about the geometric content of tensors G, 
C, n, n and their mutual relationships. What concerns kinematical content of 
formulas (|195p (|196p . compare them with the statement ([55)) valid for the gen- 
eral, non-constrained continua. Afhne framework sheds some light on geometry 
hidden behind analytical formulas of continuum mechanics. 

When the system is potential, e.g., when one deals with the hyperelastic 
afHne body, then equations of internal motion are given by (IllSp or its alter- 
native forms like pT7|) . with N'l , N^^ given by p^O)) pM|) 
(|194p . Obviously, one can admit also more general forms of the dependence of 
N, N on the configuration cp, e.g., ones describing the elastic but not necessarily 
hyperelastic affine dynamics. 

One can reasonably expect that the most useful and realistic models are 
those combining some potential term with some purely dissipative, viscous one. 



N^^ = N[VY^ + iV^,, = -I- iV^i,. (197) 



dV 

For the isotropic internal viscous friction we have 



(198) 



The last formula is written in a somewhat pretentious, but geometrically correct 
way. Obviously, n is the spatial dimension, physically n = 3, in some problems 
n — 2, — Vl\g^^ and Volo denotes the reference volume of the body. 

Traditional symbols 77, C, are used for coefficients of linear and isotropic 
internal friction. The square-root-term reduces to unitary when orthonormal 
bases are used in V, U . Geometrically the square root of det[(7,;j] is the scalar 
density of weight one in V , the square root of det[?7^B] is the scalar density of 
weight one in L/, and det \_^^a\ has a double structure: it is scalar density of 
weight minus one in V and scalar density of weight one in U . The total product 
of determinant expressions 
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is a scalar quantity, as it should be. The formula (|198p is obtained from (|123p 
when the textbook formula for the linear and isotropic stress tensor of viscous 
friction is used, 

<i = 2,7 d}^ + (^C - ^ ) 5afc d-'g'' , d^' = \ {^'' + ^'') ■ (200) 

Generalizations to anisotropic and non-linear models are obvious. A^^^^^j must 
be then some anisotropic or/and nonlinear tensor function of d^^ . 

Another model of dissipative affine torque should be used in problems of 
external friction, e.g., when one discussed an affine motion of suspensions. The 
simplest models are ones linear in f2, 

= _F-iki nki, ^ki = 9km ^"^1, (201) 

where F is a constant fourth-order tensor. In isotropic models we have 

N'^ ^-kn'^ - I n^' - p n^ag'^ , (202) 

i.e., 

pijmn ^ kg^'^gjn ^_ IgJ^gi^ + pg'^ g"^"" . (203) 

Nonlinear modifications are structurally obvious. Let us notice that the metric 
tensor may be partially eliminated by putting: 

N'^ = -F%r or N', ^ ~F)„r f^™„, (204) 

i.e., in the isotropic case, 

N) = -kn)-i ni - p n\ , (205) 

FV„" = kd\nSf + lg,^g"'+pS',S^^. (206) 

Obviously, in the mentioned applications it is rather natural to expect that 
it is mainly the g-skew-symmetric part of f2 j, i.e., (|30p. that contributes to the 
external friction. Then we have 

^ j — ~^ J : ^ jm — —g gmb ja ■ [^^ ' ) 

In the isotropic case this means that in (|205l) (I206|) we have: I = —k, p — 0, 
therefore, 

N) = -k {n'j - ni) ^ -k {n'j ~ g'^^g^t n''a) = -fc uj'j, (208) 

where to is the angular velocity ([50)) . 



6 Symmetries and conservation laws 

The complete, systematic description of symmetries, conservation laws and their 
mutual relationships is based on variational principles and Hamiltonian formal- 
ism. Nevertheless, many partial results may be obtained within the more general 
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Newton-d'Alembert framework, including the dynamics of non-conservative sys- 
tems, in particular dissipative ones. Obviously, even if not used explicitly, the 
Lagrangian and symplectic concepts are always somehow hidden behind the 
treatment. 

Let us begin with translational invariance. In mechanics of afhne bodies the 
problem of material translational invariance (material homogeneity) becomes 
diffused, and as a matter of fact, it disappears. There are two reasons for that. 
The first one is that on the level of dynamics one deals with global quantities 
obtained as mean values, integral averages performed over the material space. 
The second reason is that all formulas we use, in particular the one for kinetic 
energy, are expressed in Lagrangian coordinates vanishing at the center of mass. 
This fixed point and finite size of the body break translational symmetry in the 
material space. 

Obviously, translational symmetry in the physical space is still well-defined. 
Equations of motion (|114p (jllSp are invariant under spatial translations when 
the total force and affine torque N^^ do not depend on the spatial position 
X € M. Let us observe that this does not imply the conservation of kinematical 
linear momentum even if do not depend on internal generalized coordinates 
ip\, on their generalized velocities and on the time variable t. Indeed, if they de- 
pend only on translational velocities v'^ = ^ , equations of translational motion 
are translationally-invariant in M, but as seen from (|116p . A:* is not a constant 
of motion. It becomes a conserved quantity only if the total force vanishes, 

= 0. If is constant but non- vanishing (homogeneous field of forces), 
then equations of translational motion are translationally-invariant, however 

is not a constant of motion, either. This is obvious, because the center of 
mass motion is then uniformly accelerated. Nevertheless, there is some explic- 
itly time-dependent constant of motion somehow corresponding to translational 
symmetry. This is 

- fH. (209) 

Obviously, is the initial velocity at the time instant t — 0. Galilean boosts 
do not preserve equations of motion, but they preserve their general solution, 
i.e., the set of uniformly accelerated motions (and separately preserve the subset 
of uniform motions). There is some prescription which associates with this 
symmetry some time-dependent constants of motion, namely: 

iJ pi 

C^x'--t+ — t'. (210) 
m 2m 

It is seen that those are initial coordinates at the initial time instant t — Q. The 
simple formulas have a very interesting interpretation in terms of symplectic 
geometry and Hamiltonian mechanics, however, there is no place here for a 
more detailed discussion. 

Much more interesting are problems concerning dynamical symmetries of 
internal degrees of freedom. To simplify discussion we neglect translational 
degrees of freedom and considerate internal dynamics as autonomous one. 

Let us now begin with the spatial internal transformations, i.e., with (I46p 
where we put: b = Idjj. Equations of internal motion, i.e., IjllSp and their 
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byproducts, are invariant under such transformations, i.e., their general solution 
is transformed onto itself if and only if the affine torque satisfies the following 
transformation rule: 



(211) 



i.e., the prescription for iV as a function of state variables is "transparent" under 
the left action of any a E GL{V). In terms of the co-moving description this 
simply means that 



^-^„.,„f.«)=».«(..f,<J. (212, 

If such a rule is to be satisfied for all a e GL{V), then, obviously, N must 
be built exclusively of quantities with the capital (material) indices. It is clear 
that without additional geometric objects, such a prescription does not exist. 
The only purely material objects we have then at disposal is 51 and its tensorial 
byproducts fi™, m being non-negative integer; there are also invariantly defined 

scalars Tr (^^^^ • But the only second-order tensors built of D, are mixed ones 

in U, whereas N must be twice contravariant. The second index of N"^b may 
be obviously raised with the help of 77, the material metric, 

N^^^N^cV^B; (213) 

but one must be aware that the occurrence of 77 restricts the material GL{U)- 
symmetry to isometrics 0{U,rf). But well, we are fighting now for GL{V)- 
symmetries. Therefore, any material tensors might be formally fixed and used 
for producing N^^ from VL^Li e.g., in the simplest case we might use the scheme 



N 



AB 



j^AB^L l^y^^P^ flK^^ ^214) 



etc. Obviously, such strange models of N are completely useless for describing 
the elastic-like behavior. For such purposes we need the Green deformation 
tensor Gkl ■ Incidentally, let us remind that G may be also used for shifting the 
material indices, like r] in (I213P although this is not always physically motivated. 

Summarizing: Realistic dynamical models compatible with the idea of in- 
variance under spatial transformations have the form: 



^AB 



^ N^"" (Gkl^^'^d) . (215) 



Incidentally, let us remind the isotropic constitutive laws for unconstrained 
continua, when the second Piola-Kirchhoff stress tensor is expressed as a func- 
tion of the Green deformation tensor and the material representation of the ve- 
locity gradient. Their analogy to (I215P is obvious and certainly non-accidental, 
namely, in virtue of (|123[) . (|215p is obtained as the material average of the 
unconstrained constitutive law. 
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The spatial invariance is then automaticahy reduced to g-isometries E{M, g) 
because G is algebraically built of the spatial metric tensor cf. ([66]) and the 
preceding comments. In the hyperelastic case, when the formula (|193p holds, 
the potential energy of V of the internal 0(1/, 5)-invariant dynamics is given by 
some function of G, F = W(G). Then, roughly speaking, N is the derivative of 
W with respect to G, 

N{^) - ~2DgW, N^"" = - J?^, (216) 

or, more precisely {Gab are not independent variables because G is symmetric), 



A>v(G + ..) 



= \N^''eAB. (217) 



for any symmetric e. By the very construction, ISf^^ obtained in this way is 
symmetric, and, obviously, so is N^^ . Therefore, (|117p implies that spin (internal 
angular momentum) is a conserved quantity, 

|.l(K«_A-*.).0, (218) 

If N is derived from the potential V , like in (|214p . this conservation law is 
a consequence of the Noether theorem. Indeed, equations of motion are then 
derived from the Lagrangian 

L..,.^T,.,-V„,^\j„i§i^-WiGy (219) 

And this Lagrangian is invariant under all transformations ip aLp, a G 
0{V,g); the resulting Euler-Lagrange equations are so as well, however the in- 
variance of Lagrangian itself is something more. And it is just the invariance of 
Lagrangian that implies the conservation of angular momentum. The invariance 
of equations of motion alone, i.e., condition (I215|) does not imply spin conserva- 
tion. Nevertheless, if N^^ is symmetric, this conservation law is satisfied, even 
if N is structurally non- variational, even if dissipative forces occur. It is seen 
that the relationship between symmetries and constants of motion is a rather 
delicate matter when beyond the variational framework. 

It is seen that Noether theorem excludes higher spatial symmetries than 
the isometry group. This is because the metric tensor gij is explicitly present 
in Lagrangian, both in the kinetic energy and in the potential term, where it 
enters via the Green tensor Gab- 

Obviously, to construct any explicit prescription for the dependence of N on 
G and fl in (|215p . one must use some constitutive tensors in the material space 
U. As a rule, this restricts the a priori material group GL(U) to some proper 
subgroup. 

A typical example is the anisotropic nonlinear Hooke law, 

j^AB ^ C^BKLj^^^^ ^ ]_(jABKL (g^^^ „ ^ (220) 
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It is linear in E, if C is configuration- independent, but it is explicitly non- 
linear in generalized coordinates (^V- There is an obvious analogy with the 
constitutive laws for anisotropic continua subject to large elastic deformations, 
e.g., polymer media. However, usually in such realistic nonlinear models one 
prefers rather the situation when nonlinearity a appears already on the level of 
the very relationship between E and N. 

The relationship (|220p between E and N is invariant under such material 
mappings which preserve the tensors C, 77. They form a subgroup of the material 
orthogonal group 0(C/, 77). 

The simplest possibility is, however to rely only on the material metric 77. 
Then (|220p becomes the isotropic "nonlinear Hooke law", and C is given by: 

C^BKL ^ _^ ^r,^B^KL_ (221) 

Such a model is invariant under 0{U, rj) and it is the simplest model with 
the symmetry group 0{V,g) x 0{U,jj) acting through (|46| . 

It was mentioned in this paper that in certain formulas the material indices 
are shifted with the help of the Green tensor G[tp]. One can modify (|220l) 
along such lines, and namely replacing the constant (configuration-independent) 
constitutive tensor C by the following configuration-dependent tensor C: 

^ABKL ^ ^q-1AKq-1BL ^ ^Q'^AB q-IKL _ (222) 

We do not quote the corresponding formula for N . In spite of the formal simi- 
larity of (|22ip , p22p , the second model is completely different from the previous 
one, in particular, its nonlinearity is much stronger. 

It is convenient to use the mixed tensor G (j20p . Being linear mapping in U , 
the tensors G, may be multiplied by each other and give rise to monomials like 
G'^n^GTL'^, etc., where a, 6, c,d, etc. are integers, non-negative ones when used 
with fi. The traces of those monomials are scalars in [/, some 0(C/, 77)-invariants. 
Combining the monomials with coefficients depending on the mentioned scalars, 
one obtains some mixed tensors, elements oi U ®U* = L{U), correctly defined 
as functions of G, SI. Then raising their second indices with the help of rj^^ 
or G^^, one obtains twice contravariant material tensors, elements oi U ®U , 
just the iV"*^(G, ri) in (|215p . Seemingly one might think about infinite series 
defining (|215|) . however, the number of essential monomials will be finite in 
virtue of the Cayley-Hamilton theorem. 

The class of models described above and based on using merely the material 
metric rj for producing N from G and 17 in (|215p . is distinguished among other 
ones by the very geometry of the material space. Nevertheless, other constitutive 
material tensors are also admissible. The peculiarity of the situation when the 
prescription (|215p uses only the material metric 77 as a tool for "gluing" G, Vl 
into N is that one deals then simultaneously with the spatial g-isotropy and 
material 7y-isotropy. 

Let us now just consider the problem of material affine symmetry and cor- 
responding conservation laws. As mentioned, the only thing to be discussed is 
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the action of center-afEne material group isoraorphic with GL{U). Again, if we 
start from the invariance of equations of motion, without any reference to La- 
grangian (existing or not), we conclude that the general solution is transformed 
onto itself if and only if the following two conditions hold: 



b'^Mb'^Nr"'' = J"^, (223) 
N-(^b,'^b,t] = N^^(^,'^,t], (224) 



dt ' J V dt 

for any b G GL{U). This is the invariance under (pS)) with a = Idv- Therefore, 
b must be J^^-orthogonal, b e 0{U, J~^), and N must be a function of J[<^], 
VL, where JYp] is given by ((90|) and represents the internal inertia with respect 
to the space-fixed reference frame, 

N'i ^ N'^ [Jlip]"^ . (225) 

In the hyperelastic case this means that the internal potential energy depends 
on if through JYp]] V — V(J[(^]). Obviously, the only possibility to construct 
scalars from the twice contravariant tensor J[(^]"^ is to use some other tensor in 
V, e.g., twice covariant one (the simplest possibility). Similarly, in the elastic 
case one must use some additional tensor objects in V to construct N^^ from 
JlipY^ . Obviously, the most natural possibility is just g itself, the metric tensor 
of the physical space. 

Usually we are interested in problems of spatial and material isotropy, i.e., 
invariance of equations of motion under (|46|) with a, b running over the group 
0(V,(7), 0{U,r]) respectively. The material isotropy of equations of motion is 
possible only when the inertial tensor is isotropic, i.e., 

jAB^Jj^AB^ (226) 

Then, obviously, in (|225p J[(p] is to be replaced by the inverse Cauchy tensor 
C[(^]^^, thus, we can write: 

N^^ ^N'l {C[^]^b,n%). (227) 
In the case of hyperelastic body we have, in analogy to pi6|) : 

dW 

N{^) = -2DcV, N'^ = (228) 

the potential energy V depends on ip through C, V{(p) = W(C). 

Let us summarize the above invariance analysis of affine dynamics: 

1. Equations of motion (jll5p are invariant under internal spatial rotations 
when the affine torque N is an algebraic function of G, Q, 

N = N (ch) . (229) 
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More precisely, one should explicitly insert into this expression some con- 
stitutive material tensors Cjj in the space U, 



{G,h;Cu). (230) 



Those tensors are algebraically necessary for producing the quantity N 
from ^G, . And physically they give an account of the structure of 
internal interactions. 

2. Equations of internal motion (|115p are invariant under material rotations 
in U when the inertial tensor J is isotropic, i.e., (|226p holds, J"^^ = Ir]^^ , 
and the Eulerian torque is algebraically built of C, 

N^N{C,n). (231) 

And again some spatial constitutive tensors Cy in V are algebraically 
necessary for prescribing A^ as a function of C, il, so more precisely, we 
should write: 

iV = A(C,0;Cy). (232) 

3. Equations of motion are simultaneously spatially and materially isotropic 
when both (|23ip (|232p hold and are equivalent. But this means that in 
(|230p Cy is built algebraically of 77 and in (|232p Cy is built algebraically 
of g, thus we have two equivalent representations: 

N ^ NiG,h;7j), N = NiC,n;g), (233) 

or, alternatively, 

N^N{G.n;v), N^N{C,n;g), (234) 

where, as usual, G, G denote the mixed tensors obtained from G, G re- 
spectively by the ry-shift and g-shift of indices. Unfortunately, there is 
some disorder in notation, because characters are missing. For example, 
the "roof symbol is used to denote the co-moving representation, but at 
the same time, it is also used for the metric-based shift of indices. It 
would be perhaps better and certainly non-ambiguous to use the symbols 
like ^G, for the 77-raising and g-raising of the first index, 

(VG)^B V^^'GcB, i'Cy, :- g^'^Ckj, (235) 

and similarly A'^, A^g in the same sense, but also, e.g., E'', S^, il'', fl^ , 
etc., for the metrical raising of the second index, and similarly for the 
metrical lowering, 

:= S\/^ f7«:=f7*fe/^ (236) 
A/b A-^ScB, Ng^j-.^N'^'guj. 
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Obviously, the more indices, the worse nuisance and sometimes it is more 
convenient to use fewer symbols and to comment them in words or addi- 
tional formulas, cf. back to the comments to formulas (|19ip (|194p . 

But sometimes it is more convenient, or perhaps just necessary to be more 
pedantic. It is so even in spite of the often used argument that one can 
always simply use orthonormal bases in which r]AB —* Sab, gij =* ^ij, 
and the last formulas become trivial. 

Though we witnessed above the situation where, in a sense, indices are moved 
in the sense of Green or Cauchy deformation tensors. For various reasons not 
always orthonormal bases are most convenient, moreover we often must work 
in curvilinear coordinates, and also in curved manifold. But even if we use 
orthonormal bases, we can easily commit mistakes when forgetting about the 
geometric status of second-order tensors represented analytically by matrices. 
Any non-singular quadratic matrix a with real entries induces three transfor- 
mation rules in the linear spaces of real matrices of the same order: 

X ^ axa^^ , X ^ axa^ , x ^ a^^"^ xa~^ ^ (237) 

corresponding respectively to mixed tensors, twice contravariant tensors, twice 
covariant tensors. In complex algebra there are also other rules, first of all: 

X ^ axa^ ^ X ^ a~^'^ xa'"^ . (238) 

So, let us repeat more carefully what we said above about the structure of 
affine torque for the doubly invariant models, using conventions (|235p (|236p and 
also the following ones concerning the metrical transposition of mixed tensors 
X €U®U* ~ L{U), Y eV(E)V* L(y): 

{X^fB = iibcX^bv''^ = Xb^, (239) 

Obviously, those are metric-dependent operations, to be quite pedantic, we 
should have written X'^'^^\ Y'^'^^K^ however, we avoid the crowd of symbols 
and always keep metric tensors ry, g implicitly assumed. 

For velocity- independent internal forces, e.g., elastic ones, the doubly-isotropic 
affine torques N ((7-isotropic in physical space, 77-isotropic in the body) have the 
following algebraic structure: 

ji-i 

iV„ - ^;,(/Ci,...,/C„)(''G')", (240) 

a=0 

where, obviously, the matrix exponents are meant, and ICi are basic deformation 
invariants, e.g., in the form ICi — Tr (''G*); la are scalar functions of /C^. 
Analytically: 

n-l 

^/s = ^«(^) :'G^c"GS...''G% (241) 

a— r , 

a lactors 
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= V L{IC) n^^Gcun'"'' ■ . . ^'^^G^Lr/^^ . (242) 
t^o ' ' ' 

a factors G 

This resembles some known constitutive rules used in continuum mechanics, 
as expected in view of (|123p . Nevertheless, the above formulas may be derived 
without averaging continuum mechanics expressions. They apply also to discrete 
affine bodies like molecules (e.g. fuUerens) in an appropriate approximation. 
Obviously, from the naive point of view one might have expected the infinite 
series in (|240p . but of course, this series compresses to the finite sum in virtue 
of the Cayley-Hamilton theorem. Because of the same reason the summation 
may be extended over any range of integers from some to g Z to (to + n — 1). 
Obviously, this changes the functions la- As usual in matrix calculus, for any 
square matrix X, the zero-th exponent is taken to be identity matrix X° = I. 
When written in terms of Af -spatial geometric objects, the formulas (|240p (|243p 
(|246p become respectively: 

n-1 

Ng = Y.^- (^1' ■ • ■ '^") KT ' K'T^ = (243) 
or, equivalently, 

n-1 

n;,^Y.^- (^) ^s'^kC^'^^n ■ ■ . c;'^j , (244) 



i. facto 



I.e., 



-1 

N^^ = yia{IC)C-''''g,^C"'"...C-'-'g,,C~'''. (245) 



a=0 



a factors C 

If the internal forces are derivable from some potential depending only on 
deformation invariants, 

y = [/(/Ci,...,/C„), (246) 

then the above formulas are characterized by the special form of controlling 
functions la, namely, one can shaw that 

dU 

la{1C) = -2a—. (247) 

diCa 

From the formal point of view, (|243p is much more general than (|247p . be- 
cause nothing like the vanishing of "curl" 

is assumed. But, one must mention, quite often some doubts are raised if 
non-potential velocity-independent forces are physically realistic (although, of 
course, mathematically well-defined). 
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Let us go to velocity-dependent afHne torques. We mean the dependence 
on internal velocities. Everything expressed in the formulas (|229p (|230[) (|23ip 
((^5^ ((^55)) The analogues of ([MO)) are more complicated. Namely, 

(|240p is replaced by a sum of monomials like 

la [inv ("G, f]^'")) ) ("G)" (f]^'")) " . . . ("G)^ " . (249) 

In words: We take some products of linear operators CG)" 51^ 
multiply those monomials by coefficients depending on scalar invariants built 
of ''G, h, and take the sum of resulting expressions. Obviously, the 

mentioned scalars are traces of operator monomials. The exponents at are 
integers, and those at Q, iV'-^'> are non-negative integers. As usual, it is sufficient 
to take exponents from the range (0, . . . , (n — 1)); this is a consequence of the 
Cayley-Hamilton theorem. 

Similarly, p43p is replaced by a sum of the corresponding operator mono- 
mials in L{V), 

la (/nw(^'G,i),r!^(3))) i^cfn'^ (^!^(9))^..(9G)^^]^ (ri^^^))". (250) 

Literally meant expressions like (|249p (I250p are, so-to-speak, "neurotically" 
too general. They represent what algebra does answer to the inquiry concerning 
the most general doubly isotropic prescriptions for the afiine torque. But this 
is pure algebra, in physics only some special simple models are realistic. 

In applications the affine torque N, just as its co-moving representation 
very often is given as the sum of two terms: one depending only on the con- 
figuration ip, and the other one describing generalized forces which depend in 
an essential way on velocities. The first term describes in particular the elastic 
and hyperelastic behavior. When one deals only with purely internal interac- 
tions, it is (7-isotropic, when the constitution of the object does not distinguish 
any "material" direction, it is ?7-isotropic. These two situations are described 
respectively by ([^^ (P^O)) . but without the 0-dependence and by (^^1) (^5^ . 
so respectively 

N^NiG;Cu), N^N{C;Cv). (251) 
For the doubly-invariant models this reduces to: 

iV = iV(G;?7), N^N{C;g). (252) 

and this may be alternatively written as (|240p (|243p . If the system is potential, 
e.g., hyperelastic, then ([247| holds. 

What concerns velocity-dependent forces, the simplest model in the special 
case of discretized continua was p98l) , on the basis of linear isotropic viscoelas- 
ticity (|200p (|123p . This does not seem adequate when one deals with discrete 
systems, and then it is rather more natural to postulate the internal friction 
term in affine torques in the form: 

N,nt.d^.s = ~a + na^) - p Trin) g-\ (253) 
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or using the notation ([50)) (|236p . 

N^nt.d^ss g - (f7 + n^^^)^ - /3 Tr{n)l = -2ad - P Tr{d)I, (254) 

obviously, I denotes the identity operator in V, and a > 0, /? > 0. 

Let us notice that in the physical case of weakly compressible (almost iso- 
choric) objects, the formula p98l) contains higher-order terms, nonlinear in state 
variables. It is seen that (|253p / (|254p is but the very special case of (I250p . Using 
the technique of invariant tensor expressions we can write the general expression 
for the isotropic fluid-type internal dynamics in the form: 

n-l 

Nint.diss g — fa ■ ■ ■ ^n) d"" , (255) 

where Ca are scalar invariants built of d according to the standard trace pre- 
scription, 

Cb^Tr{d^). (256) 
The formula (|255l) may also contain the non-dissipative pressure term 

Nprg = -pln, i.e., Npr = -p g^^. (257) 

Obviously, if the internal friction is anisotropic, the prescription for N as 
a tensorial function of d, must contain some constitutive tensors, e.g., in the 
linear case: 

NlLd^ss^-^'^'^'d^b, (258) 

the shift of indices meant in the metrical g-sense. 

It is easy to reformulate the above expressions into language of co-moving 
geometric objects (tensors in U). 

Finally, let us mention about other, very important dissipative problems, 
namely, the external surface friction. Applications are obvious: imagine a ho- 
mogeneously deformable small suspension or inclusion moving in fluid. It is 
not only translational motion but also the internal motion in iy9-degrees of free- 
dom that is faced with frictional obstacles, we mean the friction between the 
surface of "suspension/inclusion" and the surrounding medium. As usual, the 
simplest and most natural assumption is that of generalized friction forces linear 
in generalized velocities, e.g., in the isotropic case, 

Next.dtss g = -vVt. (259) 

This simple expression, however, looks rather not very adequate, because it 
is physically natural to expect that the internal motion is obstacled different 
way in the special cases of rotational, shear-like and dilatational motion. So, it 
is reasonable to suppose something like 

Nlt.Mss = - P [d-' - ^ 9a, d^' g'^^ - 7 9ab d^'g^^ ^ 

= -auj'^ ~ /3d'^ - (7 - 9at d'^'g'^, (260) 
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where a, /3, j are positive constants, and the meaning of symbols w, d is Hke in 
pO|) ((3T|l . Obviously, one can also discuss anisotropic models when the above 
constants a, /3, 7 are replaced by some fourth-order tensors. Another a pri- 
ori possible modification is the external friction nonlinear in velocities. In the 
isotropic case this will be again obtained from the combination of operator 
monomials (rfg)", a = 0, 1, . . . , n — 1 with coefficients depending on the scalar 
invariants of d, i.e., quantities Tr{(fi), a — 0,1, ... ,n — 1. 

7 Towards afRne dynamical symmetry 

Two very important and at the same time very delicate points were stressed 
many times in this paper. Let us repeat and discuss them, to be able to finish 
this step of investigation with some some conclusions opening the perspective 
on further developments. 

1. There are two ways of deriving equations of motion of complex and con- 
strained systems: 

(a) the procedure based on Newton equations and d'Alembert principle 

(b) the procedure based on the variational Hamiltonian principle, Hamil- 
tonian formalism and Poisson brackets. 

2. Kinematic of our system is based on affine geometry, however, its dynamic 
in not invariant under the action of afhne group. The highest dynamical 
invariance we were dealing with above, was that under isometry groups, 
both in the spatial and material sense. Therefore, there is only partial 
analogy between our equations of motion and gyroscopic Euler equations. 
This is at least aesthetically non-satisfactory and disappointing, and brings 
about some questions concerning the status of dynamical afhne symmetry 
in mechanics and fundamental physics. 

The procedure [Ta] in the item [T] is rather more popular among specialists in 
continuum mechanics. It is directly applicable both to conservative and dissi- 
pative systems. The balance form of equations of motion appears there in a 
rather natural way , nevertheless, there is no direct and systematic relationship 
between symmetries and conservation laws, only some intuitive hints do exist. 
Unlike this, the procedure[Tb]offers a systematic theory of that relationship; it is 
based on Noether theorems. On the other hand, dissipative terms of equations 
of motion are then introduced "by hand", and they are always more or less ex- 
ternal, exotic corrections, a "foreign body" in Hamiltonian framework. But the 
balance form of equations of motion appears there simply in a canonical way. 

The item[5]is strongly related to that problem. Namely, the search of affinely- 
invariant dynamical models is much more easy within the Hamiltonian approach 
with its direct relationship between symmetries and conservation or balance 
laws. Namely, it is seen that the first and main obstacle against dynamical 
affine symmetry is due to the position of the spatial and material (reference) 
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metric tensors g, rj in the "usual" expression for the kinetic energy. In other 
words, it is due to the very particular Euchdean structure of the configuration 
space, namely, one implemented by Euclidean structures in M, N, the physical 
and material spaces, according to (jllOl) . this metric T is given by 

r = mg®{g(E)J) , T'^ = — g-'^ ® (g^^ (E) J'^) . (261) 

m 

However, from the purely geometric point of view other metrics on the config- 
uration space are much more natural, ones partially or completely independent 
on metrics in M, A^.This is geometry, but some physical motivation was also 
outlined in the Introduction; in any case, physical models may be formulated 
only when mathematical background is prepared. 

The apparently strange formulas (I176P (|177p . i.e., an alternative expression 
of (|110p is a good starting point. In (|110p we were dealing with a quadratic 
form of generalized velocities with constant coefficients built of g and J. Unlike 
this, in (I176p (jl77|) kinetic energy is expressed as a quadratic form of geometri- 
cally nicely-interpretable non-holonomic velocities, however with configuration 
-dependent coefficients built of {G[(p],J) or {g,J[ip]). So, the bad alterna- 
tive: either constant coefficients but representation of velocities non-adapted to 
geometry of the problem, or conversely-geometric affine velocities but variable 
coefficients. Why not to take the "good" features of both schemes and just to 
unify, join together their advantages? There are two possibilities of expressions 

quadratic in (■y^llj) (jj'^^fl'^B^ with constant coefficients. 

The first one consists in replacing Gab in (|176p (|177|) by 77: 



Ttr = '-^r^ABV^v^. (262) 
1 



T„t - -t^ab^k^'^lJ'^''. (263) 

The second possibility fixes the spatial metric g and some additional spatial 
tensor h G V* ®V*] the latter one is substituted instead the configuration- 
dependent J [(^] , 

Til 

Ttr = J9^JV'V', (264) 

- l9,^.nW^. (265) 
Let us observe that (|262p (|263p may be alternatively written as follows: 

This is like (jllOp : the difference is that the usual metric tensor g is replaced 
by the Cauchy deformation tensor C[Lp]. Because of this there is something like 
the mentioned similarity to the concept of effective mass in solid state physics. 
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where: 



Similarly, (|264p (|265p may be written down as follows: 

T.„, = -5,,— ^M^] • (269) 

h[ipf^ ^ ^-^^,ip-^^j h'^ (270) 

is the co-moving, thus, configuration-dependent representation of the spatial 
tensor h. 

Particularly interesting is the special case of maximal symmetry under isome- 
tries in U, V respectively, when 

J^^=l'n^^, h'^^lg'^. (271) 

Then we respectively obtain for (|263l) (|265l ) 

T.nt = \vab^\'^''lii'^''. (272) 
I 



^g,jn\mg^'. (273) 

An important property of the model (|262p (I263P is that it is invariant under 
the total afhne group GAf (M) in the physical space. What concerns material 
invariance, (j262p is invariant under orthogonal group 0{U^r]), and (j263p is in- 
variant under 0{U,r]) n 0{U, J~^). The latter group becomes simply 0{U,r]), 
when then internal tensor J is isotropic, i.e., the first of equations (|27ip holds. 
As seen from (|266l) (|267p . the corresponding metric tensor F of the configuration 
space is given by 

r = m c[v3] © (c [(^] ® J) , r'^^ — c[ip]-^®(c[ipr^®j-^). (274) 

TO V / 

Unlike (|26ip . it is curved; the corresponding geometry in the configuration space 
Q is essentially Riemannian. It has a large isometry which contains GAf{M) 
acting through [34] with B = Mn, and 0{U,ri) n 0(C/, J~^) (in particular just 
0{U,ri) when J — It]~^) acting through (|46ll on Qint and trivially on M. 

Quite symmetrically, the model (|264p (|265p is invariant under the total 
GL{U) acting through P5|) on our configuration space and it is also invari- 
ant under E{M,g) r\E{AI,h) (in particular, under E{M,g) when h = lg~^) 
acting through (p4| with B ~ Id„. The corresponding metric tensor F on Q is 
given by 

r = mg®{g(g)h[ip]) , T'^ = — g'^ ® ( g~^ ® h[(pr^) . (275) 

TO V / 

The internal parts of metric tensors (|274p (|275p . i.e., the corresponding met- 
ric tensors on Qint do factorize into tensor products of V- and [/-terms, just 
like dM]). 
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But if we once dare to give up the d'Alembert form (jllOp there is no reason 
any longer stick to factorization. The most general kinetic energy (Riemannian 
metric) on Q splitting into translational and internal parts, and affinely invariant 
in M (invariant under GAffM acting through (|34p with B = Id^) has the form: 

T = Ttr+ T,„t = y VAB i^v^' + l C.^'ac ^\ ^""d , (276) 

where the coefficients £ are constant; obviously they are components of some 
fourth-order tensor in U . Being coefficients of a quadratic from of they are 
symmetric in biindices, i.e., 

L^'I'c = CPc'^A- (277) 
It is clear that the metric tensor F underlving 12761 has the form: 

V^mC {ip\^ dx' ® dx' + C^A^c <^"' ^3 d^'B ® d^^D- (278) 

It is curved, in none generalized coordinates its components may become 
constant. Let us observe that the dip dip-pait is autonomous, unlike this, the 
dx €5 dx-part is (^-dependent. 

Similarly, for kinetic energies afhnely-invariant in N we have: 

T = Ttr + T,nt = y 9^J vV + i TVL n\ , (279) 
thus, the corresponding metric tensor F on Q has the form, 

F = m gij dx"- ® dx^ + ip~^ ip^^ ^/ TZ^L dip^A ® dtp'^B- (280) 

This time the internal and translational part are mutually independent. Obvi- 
ously, 7?. is a constant fourth-order tensor in V, symmetric in biindices, just like 

There is no model of kinetic energy, i.e., no Riemannian structure on Q 
which would be affinely-invariant simultaneously in N and M. The reason is 
that the affine group is not semisimple and its translations subgroup is a normal 
divisor. For example (|276p (I278P is afhnely-invariant in M but its maximal 
group of A^-symmetries consists of isometries in U. And conversely, (|279p (|280p 
is afhnely-invariant in A^, but in M it is invariant at most under isometries. The 
mentioned situations of maximal two-side symmetry occur when the tensors £, 
TZ are algebraically built respectively of {r],Idu) or {g,Idv)- 

For kinetic energies (metric tensors on Q) affinely-invariant in M and only 
isometrically invariant in A^, we have 

£Vc = I VAC r?^^ + y S^cS'^A + f SVc. (281) 

/, A, B denoting constants-generalized scalar moments of inertia in affine mo- 
tion. 
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Similarity, for models affinely-invariant in N and only isometrically invariant 
in M we have 

T^h'k = ^ g^kg'' + ^ M + I S^S'k (282) 

with the same meaning of constants. 
Therefore, explicitly we have 

T,nt = ^Tr (n^Vl) + 4 (^') + f ^K^) ' ' ^^^^^ 
T,„, = ^^r(^!^^7) + |^r(^]2) + |Tr(^!)^ (284) 

respectively for (|28ip and (I282p . Let us stress that the transposition of 57 G 
L{U), fl G L{V) is always meant in the metrical sense, respectively of the 
metric tensor r/, g: 

- T^^^rjBDn^c^^B^. (285) 
{n'^Yj = g^'^g^ini^n/. (286) 

The two last terms in (I283P (|284p are pairwise identical, but nevertheless it 
is convenient to distinguish consequently between the symbols f2 and D, under 
the corresponding trace expressions. This is not only more "aesthetic", but also 
prevents from some mistakes; though the first terms in (|283p (|284p are different, 
and there it is just forbidden to confuse ft with il. 

It is interesting that there exist models of Tint which are simultaneously 
invariant under GL{V) and GL{U). They correspond to the vanishing value of 
I , i.e., 

T,nt - ^Tr (h^) + ^Tr (n) ' = ^Tr {n^) + ^Tr {nf . (287) 

Except the singular case n = 1 such a "kinetic energy" (metric tensor on 
Qint) is never positively definite. Nevertheless it may be physically useful and 
the negative configurations may be interpreted as an alternative description of 
elastic forces, without any use of potential energy term, just within the frame- 
work of purely geodetic models. 

The first term in p87p has the signature (in(n + 1), ^n{n ~ 1)). Obviously, 
this is the main term, and the second one is a merely correction, for A = Q 
the corresponding "metric" would be strongly degenerate. For the special case 
A = 2ri, B = —2 more "generally" for the ratio A : B — n : (—1), one obtains 
the Killing "metric" on the linear group (or, more "generally", something propor- 
tional to it). Obviously, this "metric" is degenerate and has a one-dimensional 
singularity corresponding to the dilatational normal divisor of the linear group. 
For the generic choice of {1,A,B) the metrics/kinetic energies (|283l) (|284p are 
non-degenerate. There exists an open subset of triples (I, A,B) e M for which 
these metrics are positively definite (Riemannian). With exception of the sin- 
gular dimension n = 1, for all such triples I must be non- vanishing. 
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It is instructive to notice that the total kinetic energy corresponding to (|262l) 
(|283p . i.e., for the M-affine and A^-metrical models may be written down in the 
alternative forms: 

T = !^r^ABV^v^ + lvKLV^'''n'<Mn\+ (288) 

m I f, I A B , 

— ~^ Cij w + 2 Cki C il'^r^nH 2"^-' ^ 2" ^ * ' 

Similarity, for the kinetic energy corresponding to (|264ll . (|284p . i.e., one 
affinely invariant in N and metrical in M we have: 

+ ^ n'j Q-'i + ^ n'l = (289) 

= y 9^J v'v' + \ 9ki 5™" + y + | n\ , . 

A complete description of the scheme of breaking the affine symmetry and re- 
ducing it to the metrical one is achieved when some additional metric-dependent 
terms are admitted. The corresponding expression for the kinetic energy, i.e., 
for the metric tensor on Q, has the form: 

T = \{miGAB+fn2TiAB)v^v^ + (290) 
+ \ (IiGklG*'^ + I2 riKL r,''" + I3 Gkl ?7*'^ + 

+ I4 VKL G*^^) n'^M ^\ + y n\j h'l + 1 h'l h-'j 

or, alternatively, 

T = ^ ("^1 9ij + ^2 Cij) v'v^ + 

+ ^(Ii5M5"" + l2Gfc,G"")f}'=™r!'„+ (291) 

+ I (I3 9ki G'"" + 14 Cki 5'"") n'^ + 

A ■ B 

~t~ y ^ j "t" ^ i ^"^j . 

If translational degrees of freedom are active, then the two-side affine invari- 
ance is not possible. It is possible only when we formally put mi — 0, m2 — 
(translational degrees of freedom neglected), and in addition Ii = 0, I2 = 0, 
I3 = 0, I4 = 0. Then the metric tensor on Qint is affinely-invariant both on 
the left (in space) and on the right (in the body). The total affine invariance in 
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space is obtained when mi = 0, Ii = 0, I3 = 0, I4 = 0. The total affine invari- 
ance in the material sense corresponds to the choice: 7712 = 0, I2 = 0, I3 = 0, 
I4 = 0. For any choice of constants in (|290p (|29ip the corresponding kinetic 
energy (metric tensor on Q) is invariant under spatial and material isometrics. 
All those metrics are curved (essentially Riemannian), except the special case 
7712 = 0, Ii = 0, /2 = 0, I4 = 0, A = 0, B = 0. In this special case the metric 
(pgij) becomes flat (Euchdean) and reduces to ([TTU)) . i.e., p7S)l . PT7)1 . or, 
more precisely, to its particular case J"*^ = Ir]"^^ , i.e., the spherically symmet- 
ric top subject to homogeneous deformations. Therefore, mi = m, I3 = I, m 
denoting the usual mass of the body and I the scalar inertial moment of the 
isotropic top. 

It is both easy and instructive to write down explicitly the Riemannian 
metrics F on the configuration space Q, underlying the above kinetic energies. 
Namely, for (PTOl \^UT} they are given by 

F = - (mi gij + 7712 Cij ) dx^ (8) dx^ + 

+ (li .g., + I2 a, 7/^^ + I3 g^J v^"" + I4 aj (292) 

+ A ip^'^^j tp^'^^t + B p)'^^i <^"^^j) dip>\®dLp^B- 

This is the family of metric tensors on Q, ordered in a hierarchic way on the 
basis of their isometry groups. As mentioned, certain choices of constants cor- 
respond to isometry groups containing GL(y) acting on the left, GL{U) acting 
on the right, and sometimes both of them if translational degrees of freedom 
are neglected. 

More general choices of controlling parameters correspond to situations when 
the isometry groups in Q are smaller and based only on isometrics in Af , N. If 
they contain transformations induced by affine isomorphisms of M or /and N , 
then certainly the metrics F are curved (essentially Riemannian). Their inter- 
esting feature is that such affine models may describe bonded elastic vibrations 
without any use of potential energy. The elastic dynamics may be encoded then 
in the very form of appropriately chosen kinetic energy, as a purely geodetic 
motion in Q. Without affine invariance this would be impossible. In particular, 
if m2 0, Ii = 0, I2 = 0, I4 = 0, A = 0, B = 0, the corresponding metrics 
on Q are flat (Euclidean) and the general solution consists of straight-line in Q, 
evidently non-bounded, non-physical behavior. 

When dealing with models admitting hypothetic affine symmetry, at least 
partial one, we cannot rely upon the d'Alembert principle in its traditional 
formulation. The only natural procedure one has at disposal then, is based on 
variational principle and Hamiltonian formalism. Dissipative forces are then 
postulated as some correction terms motivated by some phenomenological and 
intuitive guiding hints. This procedure, based on Poisson brackets and Legendre 
transformation, was described in section 5, however it was specialized there to 
non-affine models. Let us now write some explicit formulas for metrics affinely 
invariant in AI or in N. In the case of affine symmetry in space, p76|) ([278]), 
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Legendre transformation has the form: 

PA^mr^ABV^ , S^s = fi^c. (293) 

Inverting it we obtain the following formula for the kinetic Hamiltonian, i.e., 
expression of energy through canonical variables: 

r^rtr+ %nt = ^ V^"" PAPB + ^ C^B^'d ^^A^^'c, (294) 

2m 2 

where 

^Vl C^k^'d = S^D S%. (295) 
The corresponding contravariant inverse of the metric tensor F (|278p is given 

by 

r-^ = - cM-i + ^Vi, ^'a (296) 

Similarity, for the metrically afhne models (|279p (|280p Legendre transformation 
is given by 

p, = mg,jv' , E*, =7^'J^^!'fc, (297) 
and the phase-space expression for kinetic energy becomes 

T^Ttr+%nt = ^ g'' P^PJ + ^ ^''a^'^c, (298) 

Zm 2 

where 

n\''in'k'd^s%6%. (299) 

The corresponding inverse metric has the form: 

m ox^ oxJ oip^B oip^D 

For the general models p76p p79p it is rather difficult to find the explicit 
expressions for p94p p98p .i.e.. for the inverse coefficients £, TZ. It is also difficult 
for more specified models (|290p (|29ip controlled by eight scalar coefficients. 
However, it may be easily done explicitly for the very special models (|28ip / (|283p 
and for (|282p / (|284p . and the more so for the simplest models based on (|287l) . i.e., 
corresponding to vanishing /. And it is just these particular models which seem 
to be most interesting in dynamical applications and in theoretical analysis. 
In the sector of internal variables (relative motion) Legendre transformations 
for the /:-models (P5T|) / (^55)) and 7^-models ([^ / (PM)) have respectively the 
following forms: 

E^B = C\''d^''c, (301) 

S^L = l77^*SLArr!V + Af^^L + i3^^Lf2*V, (302) 

E', = TV,\n\, (303) 
E^j = I5™gj„17"„ + A17^j+B(5'jf7"\„. (304) 
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In the sector of translational variables, we have respectively 

Pa = mrjABV^, (305) 
Pi = mgijV^. (306) 

Obviously, (|30ip (I302p p05p may be as well expressed through the spatial V- 
representation, and p03l) p04p (|306p - through the material {/-representation. 
This is, however, non- natural (although sometimes useful in a sense). The con- 
stant tensors on the right-hand sides are then replaced by ip dependent ones. For 
example, in p02ll (|305p and (|304p (|306p the tensors ij, g are replaced respectively 
by the Green and Cauchy deformation tensors G , C [(p] . 

In general it is rather difficult to inverse effectively the formulas (I30ip (|303p , 
however this may be easily done for the special cases p02l) p04p . and it is just 
these special cases what is particularly interesting both from the point of view 
of geometry and applications. One can easily obtain then the explicit form of 
Hamiltonian formalism. 

The inverses of p02p (|304p may be respectively expressed as follows: 

n\ = i^^*^r;i^.S% + iE^i + i^^i£^,,, (307) 
/ A B 

= i5"".9jnS'\n + is^-i<5^S™„„ (308) 
/ A B 

where the inverse inertial constants /, A, B, are given by: 
I=j{l^-A^),A=^{A^-I^),B = -^(I + A){I + A + nB) . (309) 

The simplest situation is when the internal kinetic energy is affinely-invariant 
simultaneously on the left and on the right, i.e., when 1 = 2. Then, obviously, 

i = (infinite/), A = A, B = -—A{A + nB). (310) 
/ V / B 



If there is no i?-correction term ,B = 0, then, similarly, 

(311) 



i = (^infinite B 

In virtue of p07p (|308p . the corresponding kinetic Hamiltonians (geodetic 
Hamiltonians) are given by: 

%nt = ir;KLS%S^wr;^^^ + iE^iE^K + is^KS\ (312) 
2/ 2A 2B 

%nt = ^g^ky^'o^\9'' + + (313) 

21 2A 2B 

When / = (i = 0) then the general solution of the geodetic problem is 
given by exponentials: 

ip{t) = exp [Et) ipo = (^0 exp [tp^)^'^ Eipot) = Po exp {Etj . (314) 
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Here E is an arbitrary element of L{V) or equivalently, E is an arbitrary element 
of L{U), and (fo is an arbitrary element of LI(U, V). Roughly speaking, they 
are constants of motion, or initial conditions : 

(^(0) = ^0 , J) (0) = E^o = ^oE. (315) 

If / 7^ 0, pi4p is not any longer a general solution of the geodetic problem. 
Nevertheless, even then there exist special solutions of the type (|314p . so-called 
stationary solutions. They are special in that initial conditions are subject to 
certain restrictions; namely, if we use the representation: 

(^(t) = (poexp(M) , (316) 
for stationary solutions of pi2p ; then E is t;- normal in the sense that: 



E,E'>^ 



0, (317) 



where E^ is the 77-transpose of E, 

{e^"^Yb:=V^^VbdE^c- (318) 

Roughly speaking, E^'^ does commute with E^ . This holds in particular when 
E is 77-symmetric or ?7-antisymmetric: 

E''^ = ±E'^. (319) 

Similarly, for stationary solutions of (|313p we have: 

(^(t) ^cxp(M)^o, (320) 
where (po € LI{U, V) is arbitrary, just like in pi6p but E E L{V) is g-normal, 

[E,E3^] = 0, (321) 

where 

{ES^Yr-^f/'m.E^k. (322) 

This type of "stationary solutions" is interesting in itself, just some curious 
counterpart of stationary rotations in mechanics of anisotropic rigid body. 

But this was some kind of digression. What is maximally interesting, these 
are doubly affinely-invariant geodetic models pi2p pi3p with / = (inciden- 
tally, they are identical in both versions of the formula). Then the general 
solution is given by the matrix exponents pi5p . Of course, one can admit in 
addition to T some potentials V and consider other models, then no longer ones 
admitting exponential solutions. But it is a very curious circumstance that even 
within purely geodetic framework one can describe strongly nonlinear elastic vi- 
brations. Dynamics is not then encoded in anything like but just in the 
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kinetic energy, i.e., in the metric tensor of the configuration space. This resem- 
bles some properties of the Maupertuis principle. More precisely, this is true for 
the isochoric (incompressible) part of motion, when we consider only degrees 
of freedom ruled by the special linear groups SL{V), SL{U). General solution 
contains then an open subset of bounded motions and an open subset of non- 
bounded, escaping and collapsing solutions. Roughly speaking we are dealing 
with some dissociation threshold and bounded, non-linearly vibrating processes. 
And all this without potential, and analytically based on the properties of ex- 
ponents! The purely dilatational part of geodetic motion is non-bounded or 
collapsing, except, of course the constant solution. But this purely dilatational 
part may be stabilized by introducing some auxiliary dilatational potential in 
one dimension, some oscillator, potential well, etc. 

All this is a very important argument for investigating afhnely-invariant 
dynamical models, although from some point of view they might seem "exotic" 
But they are not more exotic than the concept of effective mass in solid state 
physics. Analytical tools of the analysis are based on the properties of the matrix 
exponential map and the polar and two-polar decomposition of the matrix 
representing the internal configuration. 

First of all, let us notice that the expressions (|312p pi3p may be written 
respectively in the following forms; 



Tint — 



213 \ J Afi 



%nt = -Tr{^^)+-iTr^r^-Tr{S') (324) 

where, let us remind, the tensors S, V denote the canonical spin and vorticity, 
given by ([70]) ([7T|l . The constants a, /?, /i are expressed by A, B, C as follows: 

a^l + A, p ^ -{1 + A){1 + A + 7iB)/B, ^={l^-A^)/l. (325) 

The formulas (|323p p24p may be also expressed as follows: 

7-. = ^^(2) + ^C(l)^ + ^ll^f (327) 

where C{k) denotes the fc-th degree Casimir quantity, and \\S\\ are the 

magnitudes of vorticity and spin, 

\\Vf - -lTr{V'), \\sf = Arris') (328) 



C(fc) ^ Tr (S'^) = Tr . (329) 

It is seen that the difference in the invariance properties between (|312p and 
pi3p . i.e., between (|323p and (|324p is reflected only by the last, i.e., third, 
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terms in (|323ll (I324p . For more general models, like (|290p (|29ip . situation is 
more complicated and the relationship between velocity-based and canonical 
models is more complicated, although it still does exist. 

Let us stress that p26p (|327p and the doubly (left- and right-) affinely in- 
variant geodetic models (I287P of internal dynamics differ by constants of motion 
proportional to \\V\^ and \\S\^ . They are, so-to-speak, half-affine, i.e., affine in 
the space and metrical in the body (I326P or conversely, metrical in the space 
and affine in the body (|327l) . They are special cases of (|290p (|30ip . It may be 
interesting to consider their combination given in Hamiltonian terms by 

7I„* = ^C(2)4-^C(l)^ + i-rf + i-||5f, (330) 

again with a, /3, /i, v being some inertial constants. 

However, for us the most important expression would be one containing only 
the Casimir a-, /3-terms, 

%nt = ^ C(2) + i- C(l)^ (331) 
or merely, the main second-order Casimir term 



T^nt = ^ C{2) = S^ E-', = S-^sE^^. (332) 
2a 2a 2a 



It is just this simplest model where one can observe the mentioned encoding 
of nonlinear elastic vibrations in negatively-determined part of the kinetic energy 
expression. To see this explicite, one should use so-called two-polar splitting for 
the mapping ip G L{U, V^)describing the internal configuration. 

Any linear isomorphism (p oi U onto V induces in linear spaces U, V two 
metric-like tensors, just the Green and Cauchy deformation tensors G[(p] £ 
U* ®U* , C[ip\ €V*<S!)V*. But the material and physical translation spaces U , 
V are endowed also with some fixed, thus also constant in time, metric tensors 
77 S U*®U* , g G V*®V* . And, as we saw, these metrics enable one to construct 
the mixed tensors G[(p\ S U ®U*, C[lp\ G V ()^V*. They have eigenvalues Aa, 
Aa~^, a = 1, . . . , n. In certain formulas it is also convenient to use the quantities 
denoted by Q", g", a = 1, . . . , n, where 



g° = exp(g'^) = VAa. (333) 

Obviously, A^-s are positive, because 77, G[ip\ and C[(^] are positively- 
definite tensors. Using the quantities as diagonal entries, we can construct 
the diagonal matrix D = diag{Q^ , . . . ,Q"'). It may be identified with some 
linear isomorphism of M" onto M" , D : M" ^ M". 

The quantities — D— (no summation convention!) are the basic stretch- 
ings of the system. They contain the information about how much the object 
is elongated/shortened in principal directions. But they do not tell us anything 
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about the orientation of those stretchings in U and V. This information is en- 
coded in orthonormal systems of eigenvectors of G and C. The eigenvectors 
satisfy equations 

GRa = exp(2(?")i?„ CLa = exp(-2g'^)L„. (334) 

The bases Ra- La are assumed to be orthonormal, 

7] {Ra,Rb) = VKlR\R\ = Sab^g {La, ifc) = g^JL\L^. (335) 

To be more precise, they are assumed to be unit vectors; their orthogonahty 
is a generic situation for the case of simple spectra of operators G, C. 
The dual bases in [/*and V* will be denoted by i?"^, i"^, thus 

(i^^i?b) = n, (L^L,) = n. (336) 

Therefore, the deformation tensors may be expressed as follows: 
GM = ^exp(2g'^[^])i?«[v5]0i?'^[^], 

a 

(337) 

C[ip\ = ^exp(-2g''M)L'^[^]®i'^[^]. 

a 

Every mapping of internal configuration if is represented by three kinds 
of objects, namely, by a pair of metrically rigid bodies with configurations 
(. . . , i?a[(/9], . . .), (. . . , La[</']j • ■ •) in U and V , and the n-tuple of fictitious mate- 
rial points with positions (. . . , g'^[(^], . . .) in the real axis R. Those subsystems 
have respectively n{n — l)/2, n{n — l)/2, n degrees of freedom. It is important 
to note that this representation is not unique. It is so even in the special case 
of Lp with non-degenerate spectra of G[if\, C[lp\, because the permutation of 
indices a in p37p does not affect the sum. When the spectrum of G[(^], C[tp\ is 
degenerate, the non-uniqueness of (j337p becomes continuous. This singularity 
of coordinate system resembles that of spherical or polar coordinates at r = 
(yi, "d completely non-determined), although in concrete details it is much more 
complicated. 

Obviously, the linear frames L = (. . . , La, ■ ■ ■), R — {■ ■ ■ , Ra, ■ ■ •) may be 
canonically identified with some isomorphisms L : M" -> V and R : R" — >■ 
U. And, obviously, their dual co-frames L — (. . . , L°-, . . .), R — {. . . , R'^, . . .) 
correspond canonically to some linear isomorphisms L~^ : V — )> M" and : 
U — !> M". If the diagonal matrix with deformation invariants Q'^ on diagonal is 
identified with a linear isomorphism D of R" onto R", then, making use of all 
the above identifications, we can simply write 

if = LDR-^; (338) 

in matrix terms this means: orthogonal times diagonal times orthogonal. 
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It is clear that the spatial and material isometrics preserving orientation in V 
and U, A € SO{V, g), B G SO{U, 77) act on the left on the L- and i?-gyroscopic 
part of (p, because the mappings 

AL, BR (339) 

imply that ip transforms under them as follows: 

Lp^AipB-^. (340) 

Their Hamiltonian generators are respectively the spin (for A) and the neg- 
ative vorticity (for B). 

Let us now consider in a more detail the structure of L- and i?-rigid bodies. 
Their "physical spaces" are respectively V and U . And their material spaces are 
simply identified with R". The group SO{n,M) acts on them on the right, 

LC, R^ RD (341) 

C,D ^ SO{n,M.). These mapping are not expressible like p40p . by an action 
on if as a whole. Nevertheless, locally they are well-defined and so are their 
Hamiltonian generators. One can define their "spatial" (in the V, J7-sense) and 
co-moving (in the R"-sense) angular velocities and canonical spins. 

The "co-moving" and "spatial" components for the angular velocity of the 
L-top are defined as: 

L%^V„ (342) 

or x = X%La®LK (343) 

Obviously, x G SO{V,gy, x G SO{n,Ry. The first latin characters a, b are 
simply labels, i.e., indices in M", whereas the middle ones i, j are tensor indices 
in V. Similarly, the "co-moving" and "JJ-spatial" components of the angular 
velocity of the i?-top are given by analogous formulas: 

n := (r\j^R,^=R'^mj^R\, 

(^j^R.\^R%i, or ^^d^Ra^R^. (344) 

Again for this rigid body U plays the role of the "physical" space in spite of its 
being originally "material". The "material" space for the i?-top is again M" . The 
choice of non-holonomic velocities as (^q'^ ,x%,'{^%^ or (9", X j 1 ^'^m) depends on 
particular purposes. Just as is was the case with the usual, physical rigid bod- 
ies, it is convenient to introduce non-holonomic, just Poisson-non-commuting 
canonical momenta {pa, p°'b,T%) or {^pa, p^j ,t^b) ■ Here pa are canonical mo- 
menta conjugate to deformation invariants q°'. The quantities p%, f% and pj. 
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T^B are canonically conjugate to x°bi X j"> '&^m- Obviously, p,T G SO{n, M)', 
p £ SO{V,gy, T G SO{U,r])' . The meaning of the mentioned "conjugacy" is en- 
coded in the following duality formulas: 

((p,r,p),(x,^9,g)) - ((/5,f,p),(x,^,g)) = (345) 

= Pacf + \ Tr{px) + \ TT{Td) = Pa r + ^ Tr{px) + \ Trird). 

The Hamiltonian interpretation of p and — t as generators of the left- and 
right-orthogonal translations of ip (l^-spatial and {/-material ones) implies that 
they are simply identical with spin and vorticity. They are respectively (/-skew- 
symmetric and 77-skew-symmetric parts of S and S. Unlike p and r, their 
K"-material representatives p, f , generating p4ip fail to be constants of motion 
for geodetic systems and for Lagrangian models with potentials depending only 
on deformation invariants g". It is also clear that the constants of motion p, r 
are related to /), f as follows: 

p = p%La®L\ T^T%Ra®R\ (346) 

Let us stress that when dealing with the traditional, orthogonally invariant 
kinetic energy (jllOp (j26ip . the traditional form of deformation invariants 
(|333p is more convenient than which are particularly suited to models with 
affine symmetry. 

It was mentioned that the two-polar expansion (j338p is not unique. At the 
same time, it is well-known that the usual polar decomposition is unique. It is 
convenient to mention here about the polar splitting in this description. Namely, 
La [f] , Ra [f] are two orthogonal bases, respectively in the (V, g)-senBe and ([/, 77)- 
sense. Therefore, there exists the unique isometry U[ip\ which establishes a 
one-to-one relationship between them; 

LM^U[ip]Ra[p]; (347) 

U[ip\ E 0{U, rj; y, g). One can show that 

^ ^ U[^]Ay] = B[^]U[^], (348) 

where the linear mapping A[(p] S GL{U), B[(p] G GL{V) are respectively 77- 
symmetric, (/-symmetric and positive. Therefore, 

r]{x,A[(p]y) = T]{A[ip]x,y), g{w, B[ip]z) ^ g{B[Lp]w, z), 
T]{x,A[(p]x) > 0, g{w,B[ip]w) > 0, (349) 

for any vectors x, y, w, z. In the last two inequalities the vectors x, w are non- 
vanishing. 

The symmetric elements of the left and right polar decomposition are related 
to each other as follows: 

B[^] = U[v]A[^]U[v]-\ (350) 
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i.e, by the J7[(/j] -similarity transformation. In the usual matrix representation 
l|348l) has the form 

ip^UA^ BU. (351) 

The symmetric terms, e.g. A, may be diagonalized by the orthogonal simi- 
larity: 

A^RDR-\ i? - orthogonal. (352) 

Then p5ip becomes 

if = URDR-^ ^ LDR-^, U^LR-^. (353) 

Although L, R separately taken are not unique, U is so, and therefore. A, B are 
unique. 

Let us define now some objects instead p, f, which enable one to perform a 
partial diagonalization of the affinely-invariant kinetic energy of internal motion, 

M ■=-p-T, N = p-T. (354) 

After this substitution the second-order Casimir invariant, 

C(2) = Tr {f?^ , (355) 

becomes: 

C(2) - V ^ ^ + - V - 1 V f356) 

G (2) - 2^ + 2. ^2 16 ^„3j^2 £^ ■ i^56) 

a a,b ^^^^^^ 2 a, 6 ^"''^^ 2 

The first term in p56p admits a nice representation as a sum of the relative and 
over-all contributions. 



?1 

n 2n 



where p is the first-degree Casimir invariant, 

p = (7(1) = rr(E) = Tr (e) . (358) 

It is a quantity canonically conjugate to the center of mass of logarithmic 
deformation invariants. 



q 

n 



-Vg^ (359) 

T> ^ ^ 



The corresponding kinetic energy based on the second-order Casimir invari- 
ant C(2) as the main term, 

r=^C(2) (360) 



has, as seen form ([356]) a very peculiar structure. It consists of the first term 
which formally has a structure of the kinetic energy of the n-particle g"-objects 
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on M, and of the next two terms. The first of these terms, i.e., sinh ^-^-^ — 
term describes centrifugal repulsion between invariants, but the second one, i.e., 
cosh~^ ^ ^'^ , due to its minus-sign, represents a curious "centrifugal attraction". 
The strengths coefficients {M%)^ ^ are positive, so really one has to do 

with repulsion and attraction of deformation invariants. Repulsion is singular 
at the coincidence q"" = , whereas attraction is then finite. But at large 
|g° — g*" I -distances, attraction prevails if \Nab\ > \Mab\- This is the typical 
shape of "intermolecular" forces between deformation invariants interpreted as 
indistinguishable material points. Admitting the C(l), IIS"!!^, ||l^||^-terms, we 
do not change this interpretation. Before commenting it more, we mention only 
about some related models. 

It is seen that p56p p60p -models resemble the hyperbolic version of the 
Sutherland lattice. There are both similarities and differences. The main dif- 
ference is just the occurrence of attractive terms, the negative contributions 
to p56p . Their interpretation is slightly obscured because the coupling ampli- 
tudes M%, N% are not constant, they satisfy together with q°- s closed sys- 
tem of equations of motion. These equations may be easily expressed in terms 
of Poisson brackets between basic state variables {q°',pi,,M%,N%); obviously 
{r7°,Pfc} = n, {q\M'4 = {9",A^"4 = 0, and {pa,M^d} = {Pa,N'd} = 0. 
The brackets for M%,N'^d niay be easily obtained from their definition p54p 
and from the standard Poisson rules for p%, as Hamiltonian generators 
of the action of SO{n,R) x SO{n,R) through ([MT|l . The family of functions 
{q'^ ,Pb, M'^b, N%) forms a Poisson algebra. In the case of non-geodetic models 
with doubly isotropic potentials, depending only on deformation invariants g", 
equations of motion may be all obtained in terms of those Poisson brackets as 

^ = {F,H}; (361) 

for F one has to substitute q'^,pb, M%, N%. 

Obviously, this gives us only a partial description of motion, without the 
time dependence of gyroscopic variables La,Ra- But in practical applications 
like, e.g., vibrations of molecules, this time dependence is less important. In 
any case, having solved p6ip . we know the time dependence of p, f. Then, 
inverting the Legendre transformation, one can express the time dependence of 
gyroscopic angular velocities p, f. And then, finally, for the time dependence of 
La, Ra one obtains first-order differential equations for the n-legs La, Ra- 

- L\x\{t), = R\^'a{t). (362) 

at dt 

So, after solving the system p6ip for {q"' ,pb, M%, N%), i.e., after obtaining 
the main deformation description, we can in principle solve the next, less im- 
portant step and determine the time dependence of the main axes of inertia. It 
is so for any doubly-isotropic model. 

The very idea is that in (|356p (|360p . The "plus" and "minus" terms act in 
opposite directions, representing, respectively, the repulsion and attraction of 
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deformation invariants even in the completely geodetic, potential-free models. 
The question is, however, if the idea is correct, because the repulsion and at- 
traction strengths are not constants. There is, however, some special case where 
they are constants of motion. It is in two-dimensional models, when n — 2. In 
fact, in a consequence of the fact that the orthogonal group SO{2,M) is Abelian 
(because it is one-dimensional), in this special case we have: 



P b = p b, 



T b — T b, 



(363) 



and because of this, M^, N% are constants of motion (for the doubly- isotropic 
models), just as = 3%, t% = —V^ are. Because of this, the effective 
strengths {M%)'^, {N%)^ in (|356p (|360p are constant in time. The matrices 
L,R, D become 



L = 
D = 



cos a 
sin a 







— sm a 
cos a 



exp q 




cos ^ — sin ^ 
sin j3 cos (3 





cxp 



It is convenient to represent GL+{2,R) as M+5i(2,i 
"center of mass" and the displacement between q^ and q^, 



1 f 1 

and their conjugate momenta 



1 



Px 



{P2 -Pl) ■ 



(364) 

by introducing the 
(365) 

(366) 



The angular velocities x, i) have the obvious form: 

dp 



da 



-1 

1 



X = X 

The quantities M, N become 
M 

where, obviously. 



1 

-1 



N 



dt 



-1 

1 



1 

-1 



m^Pp-Pa, n^pp+pa 



(367) 



(368) 



(369) 



Therefore, for _B = 0, the doubly affine-invariant kinetic energy of internal 
degrees of freedom becomes 



T = ^ (Pl +P1 j + ^ — jT^Hi " 161 ^^iTZ 



(370) 
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Using the new variables f|365p (|366p (|369p . one obtains: 



9 ? 

P , Pa; 



7-= ^ 1^1 (371) 

4yl A IGAsinh^f 16^cosh^f' 

For a bit more general models (j^M)) ((5^ one obtains essentially 

the same, up to certain modification of constants and correction terms propor- 
tional to \\V\\, \\S\\. The essence of the problem remains the same, because 
are constants of motion. 

It is seen that in the geodetic problems like p7ip . the shear-part, i.e., 
51/(2, K)-component of dynamics admits a continuous family of bounded elastic 
vibrations. This is the part characterized by |7i| > \m\, where at large \x\- 
distances the "centrifugal attraction" prevails the centrifugal repulsion. This 
is in spite of the non-compactness of the x-variable and of S'L(2,IR) at all. 
The family of bounded motions has the dimension six, i.e., twice the dimen- 
sion of SL{2,M.). Above the threshold, i.e., for |m| > |n|, one deals with the 
six-dimensional family of unbounded motion. This is the typical threshold be- 
havior, which does exist on the purely geodetic level due to the curved geometry 
of S'L(2,M). In this way, the isochoric (incompressible) elastic motions may be 
described in purely geodetic terms, without any use of the potential. The only 
failure of this threshold behaviors on the level of the total GL(2,M) are just 
dilatations. According to (|37ip . the variable p occurs in a quadratic way with 
constant coefhcient, and its conjugate configuration variable q does not occur 
at all. Therefore, except the solutions q — const, p ^ 0, all solutions are sin- 
gular wit respect to q. They are either infinitely expanding, or contracting to 
the singularity q = 0. This is a consequence of the fact that G'L+(2,M) is not 
semisimple, being isomorphic to M"*" 51/(2, R) = exp(K)S'L(2, M). Of course, this 
is not a dangerous failure for our philosophy. The g-motion may be stabilized 
by introducing any bounding potential, e.g., oscillator, thin potential well, etc. 

Exactly the same holds for any n > 2. The difference is only that for n > 2, 
the spin quantities M%, N% fail to be constants of motion, they also vibrate. 
But the SL{2, M)-threshold behavior is a consequence of the commutation rules 
is S'L(2,M), GL(2,]R). More or less, the same holds also for other models with 
the partial affine invariance of kinetic energy, like (|283p (|284p . i.e., (|326p (|327p . 

Let us mention also about some other models more or less related to the 
above ones. First of all, let us remind the traditional d'Alembert model (IllOp 
(|26ip with the materially isotropic co-moving inertial tensor J^^ = Irj^^ . Then 
the kinetic energy is given in the two-polar representation by 

21 4^ + 81 2^ (Q. - Q^f + 81 2^ (Qa + Qbf • (372) 

Let us notice that now it is not logarithmic invariant but just that 
is convenient. Obviously, this expression contains only centrifugal repulsion 
of invariants. Therefore, p72p is completely useless as a model of nonlinear 
elastic vibrations. To obtain a viable model, one has to modify it by adding 
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some potential energy. In the doubly isotropic models, it will be a function 
93 (Q^,--- ,Q") of deformation invariants alone. Obviously, the first two 
terms of p72p resemble the Calogero-Moser integrable lattice and indeed there 
is some relationship here. 

It was mentioned also that (|370p (|356p (|360p resemble the hyperbolic version 
of the Sutherland lattice. It turns out that some relationship with the usual, 
i.e., trigonometric Sutherland lattice does exist as well. To obtain it. one should 
compactify the deformation invariants by the formal substitution of 

Q° = exp(ig°) (373) 

to the positively definite kinetic energy 

T^nt = -^Tr {n^) ^^Tr {n+n) , A > (374) 

Then GL{2, R) is formally replaced by U (n)- the compact real form of GL{n, C). 
Obviously, 

17 ==¥^-1^, ifeU{n). (375) 
After this formal substitution one obtains for Tint 

hut - 2a2^^- ^ 32A ^ sin^ 31^ ^ 32A ^ cos^ Sl^ ' ^ ' 

a a,b ^^^^ 2 a, 6 ^'^^ 2 

The relationship with the usual Sutherland lattice is obvious. Let us observe 
that now the deformation invariants run over the circle. Because of this there 
is no problem with the positive definiteness of the terms in (j376p . Because 
in the compact space with the circular topology it is practically impossible to 
distinguish between repulsion and attraction. 

There are also other important problems concerning afhnely rigid body and 
the problem of dynamical afEne invariance. They concern mainly two fields. 
One of them is the theory of systems of afhne bodies. And of course another 
one is quantum mechanics of such bodies and their systems. There is no place in 
this paper for studying them. Some primary discussion was delivered in papers 
quoted in the references below, some other aspects are to be investigated in 
future. 
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